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PREFACE 


The  AGARD/SMP  Subcommittee  on  “Aeroelasticity”  heard  four  technical  papers  at 
the  Fall  82  meeting  in  Toronto,  Canada.  Three  of  the  papers  represent  a  cross  section  of 
recent  activities  in  aeroelasticity,  covering  subsonic  flutter-clearance  procedures  in  Canada, 
transonic  flutter  research  in  Germany,  and  transomc  unsteady  aerodynamic  measurements 
in  the  United  States.  A  fourth  paper  by  Dr  Erwin  Johnson  of  the  United  States  discussed 
very  promising  results  Northrop  has  achieved  in  developing  an  adaptive  flutter-suppression 
system.  Wind  tunnel  tests  on  an  aeroelastic  model  of  a  wing  with  external  stores 
demonstrated  rapid  sensing  of  a  suddenly  violent  flutter  mode,  computation  of  desired 
control  system  gains  and  phases,  and  resulting  suppression  of  the  flutter-instability.  AH  four 
papers  indicate  the  great  strides  being  made  in  aeroelasticity  by  the  NATO  countries  and 
promise  even  greater  understanding  and  progress  within  the  next  few  years. 

Cj/^Ux - 

JAMES  J .OLSEN 

Chairman,  Subcommittee  on  Aeroelasticity 
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ABSTRACT 

The  Canadian  capability  in  flutter  clearance  oi  military  aircraft  carrying  underwing  stores  is  reviewed.  The  flight  test  facilitie 
and  procedures,  on-line  analogue  and  post-flight  digital  date  analysis,  and  analytical  flutter  model  are  described.  Some  results  for  the 
LAU-5003/A  rocket  launchers  carrying  C14  rocKets  armed  with  Mk  I  warheads  are  presented.  Frequencies  and  damping  values  ob¬ 
tained  from  strip  derivatives  and  doublet  lattice  aerodynamics  methods  used  in  the  flutter  computational  cooe  are  discussed. 
Experimental  result  from  on-line  and  post-flight  analyses  are  compared  for  one  aircraft/store  configuration. 

1.0  INTRODUCTION 

Most  certification  programs  for  Canadian  aircraft  with  external  stores  were  in  the  past  performed  by  the  manufacturer  of 
the  aircraft,  USAF  or  the  Netherlands.  Reliance  on  a  foreign  country*'  c:  arance  program  whose  requirements  were  no*  identical  with 
the  Canadian  requirements  could  result  in  Canada  having  to  compromise  some  op  its  operational  options  for  usage  of  the  program.  In 
1972  (Ref.  1)  a  study  for  a?,  "in  house”  certification  program  to  meet  Canadian  r  eeds  was  completed  and  it  was  concluded  in  the 
report  that  an  independent  Canadian  capability  for  external  stores  was  technically  feasible  and  economically  justifiable.  Such  a  capa¬ 
bility  would  free  the  Canadian  Forces  from  the  constraints  of  foreign  bias  and  allow  for  a  flexible  choice  in  the  selection  of  a  foreign- 
domestic  mix  of  aircraft  and  stores. 

The  first  of  an  "all  Canadian*’  effort  in  stores  clearance  on  the  CF-5  aircraft  was  in  1981  at  Cold  Lake  CFB.  Prior  to  that, 
clearance  on  the  CF-5  was  contracted  to  National  Aerospace  Laboratory  (NLR)  of  the  Netherlands.  The  flight  test  techniques  adopted 
by  the  Canadian  Forces  were  similar  to  those  which  had  proved  reasonably  successful  in  previous  Netherlands  tests  on  their  NF-5 
airciaft  and  in  a  CF/NLR  program  on  the  CF-5  (Ref.  2). 

This  paper  gives  an  overview  of  the  flutter  phase  of  the  stores  clearance  program.  The  major  participants  in  this  multi-agency 
effort  are  the  Aerospace  Engineering  Test  Establishment  of  the  Canadian  Forces,  Canadair  Limited  and  the  National  Aeronautical 
Establishment.  The  flight  test  facilities  and  procedures,  on-line  analogue  and  post-flight  digital  data  analysis,  and  analytical  flutter 
model  will  be  described.  The  discussions  will  be  limited  to  the  CF-5  aircraft  carrying  LAU-5003/A  stores  fitted  with  nose  cones. 

2.0  FLUTTER  FLIGHT  TESTS 

An  essential  part  of  a  certification  program  is  to  investigate  whether  the  aircraft  is  safe  from  flutter  for  all  stores  configura¬ 
tions  within  the  proposed  flight  regime.  Generally,  the  flutter  investigation  is  carried  out  both  experimentally  and  analytically.  The 
analytic  modelling  will  be  described  in  a  later  section.  The  flutter  calculations  include  a  large  number  of  possible  stores  configurations, 
and  some  of  the  more  critical  ones  are  usually  selected  for  flutter  flight  testing  so  that  the  estimated  damping  margins  can  be  con¬ 
firmed  by  measurements. 

2.1  Ground  Vibration  Testing 

For  the  safe  carriage  of  LAU-5003/A  stores  on  the  CF-5  aircraft,  Hatter  flight  tests  were  carried  out  at  Cold  Lake  CFB 
during  February  and  March  of  1981.  In  the  pre-flight  flutter  analysis,  ground  vibration  testing  (GVT)  had  been  performed  (Ref.  3) 
and  the  data  used  in  the  calculations  of  flutter  trends  which  formed  the  basis  for  flight  testing.  Since  no  suitable  ground  vibration 
testing  facility  was  available  in  Canada  for  stores  clearance  work,  the  data  for  the  LAU-5003A/C14  launcher/weapon  configurations 
on  the  CF-5  aircraft  was  obtained  from  a  joint  USAF/CF  test  program  at  the  Ai*  Force  Armament  Laboratory  at  Eglin  AFB  between 
September-October  of  1979. 

The  vibration  tests  were  conducted  in  the  Overhead  Soft  Suspension  System  in  which  four  cables  were  attached  to  the  jack/ 
hoist  adapters  on  the  aircraft  fuselage  as  shown  in  Figure  1.  Each  cable  was  suspended  from  an  Air  Spring  Isolator,  and  this  arrange¬ 
ment  provided  excellent  aircraft  stability  for  modal  testing  with  the  landing  gear  retracted.  Two  shakers  were  used  for  excitine  the 
aircraft  structure.  For  most  of  the  configurations,  the  shakers  were  attached  at  the  tip  tanks,  while  in  other  configurations  they  were 
attached  to  the  outboard  stores. 

A  block  diagram  of  the  data  acquisition  and  analysis  system  is  shown  in  Figure  2.  A  total  of  83  accelerometers  were  used  to 
measure  modal  response  and  two  load  cells  (one  at  each  shaker)  provided  force  input  data.  Most  of  the  accelerometers  were  installed 
on  the  right  wing  to  give  a  detailed  vibration  characteristic  of  the  wing.  Since  only  32  channels  of  data  could  be  analyzed  at  a  time, 
the  accelerometer  signals  were  separated  in  three  multiplexed  sets  with  a  common  reference  included  in  each  set.  The  data  could  be 
either  stored  on  disk  for  future  analysis  or  fed  directly  to  a  fast  fourier  analyzer.  The  modal  amplitude  and  phase  information  were 
sent  by  direct  line  to  a  CDC  6600  computer  for  additional  manipulation  and  graphic  display. 
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For  validation  purposes,  the  GVT  results  were  compared  with  thc&e  determined  by  the  Netherlands  for  the  NF-5  aircraft.  A 
comparison  of  the  resonant  frequencies  for  the  two  aircraft  are  shown  in  the  following  table: 


MODE 


NLR  NF-5GVT  EGLIN  CF-5GVT 


Wing  1st  Bending 

Sym 

3.69  Hz 

3.63  Hz 

Wmg  1st  Torsion 

Sym 

5.93  Hz 

5.87  Hz 

Wing  1st  Torsion 

A/S 

5.67  Hz 

5  74  Hz 

Wing  1st  Bending 

A/S 

6.67  Hz 

7.14  Hz 

OTBD  Store  Yaw 

A/S 

9.01  Hz 

9.53  Hz 

Centerline  Tank  Roll 

A/S 

10.03  Hz 

10.08  Hz 

The  results  are  very  close  since  the  two  aircraft  are  structurally  very  similar;  the  primary  difference  is  that  the  NF-5  is  configured  with 
a  maneuvering  leading  edge  flap  whereas  the  CF-5  has  the  standard  flap. 

2.2  Aircraft  and  Airborne  Instrumentation 

Figure  3  shows  a  CF*5  test  aircraft  fitted  with  centerline  tank  and  carrying  two  LAU-5003/A  stores  at  the  outboard  station. 

A  view  of  an  empty  rocket  launcher  with  the  nose  cone  removed  is  given  in  Figure  4.  The  launcher  may  be  loaded  with  from  one  to 
nineteen  C14  rockets,  fitted  with  either  six,  ten  or  sixteen  pound  warheads. 

The  test  aircraft  was  equipped  with  the  necessary  tett  instrumentation  for  determining  the  flight  conditions  which  included 
the  flight  altitude,  indicated  airspeed,  total  temperature,  position  of  elevator  and  aileron,  yaw/pitch/roll  rate  and  angle  of  attack.  For 
flutter  analysis  purposes,  eight  accelerometers  were  installed  on  the  aircraft  (Fig  5).  Four  of  them  were  positioned  at  the  front  and  roar 
of  the  tip  tanks  and  these  were  the  primary  ones  used  for  flutter  investigation.  Hie  remaining  four  accelerometers  were  located  at  the 
nose  of  the  stores  and  served  mainly  to  provide  data  in  the  vibration  of  the  stores.  All  accelerometers  were  of  the  piezoelectric  type. 

Figure  6  shows  the  airborne  data  system.  To  provide  data  acquisition  back-up,  all  data  transmitted  to  the  ground  station  was 
also  recorded  onboard.  The  recording  signal  conditioning  and  telemetry  transmission  equipment  were  housed  in  the  left-hand  nose 
gun  bay.  No  special  equipment  was  provided  for  the  stick  pulse  excitation.  Fc  r  simplicity,  this  was  achieved  by  rapping  the  stick 
directly,  either  laterally  or  longitudinally  as  appropriate.  The  aircraft  was  fitted  with  aileron  and  elevator  position  sensors,  and  these 
signals  were  used  for  triggering  the  start  of  analysis  of  the  Fourier  analyzer. 

2.3  Ground  Station  and  Instrumentation 

Most  of  the  flight  testing  was  performed  in  airspace  over  the  Primrose  Lake  Evaluation  Range  approximately  30  rmles  north 
of  the  Cold  Lake  CF  base.  This  range  was  equipped  with  cinetheodolites,  tracking  radar,  a  range  trials  control  center  and  an  auto 
tracking  antenna.  Data  sent  from  the  aircraft  was  transmitted  to  the  ground  station  where  data  acquisition  and  analysis  were  per¬ 
formed  (Fig.  7).  The  serial  PCM  telemetry  data  was  processed  through  decommutation  and  word  selection  systems  to  provide  parallel 
display  of  selected  channels  on  a  strip  chart  recorder.  Other  parallel  data  feeds  were  supplied  to  a  Hewlett-Packard  Fourier  Analyzer 
and  the  spectra  could  be  displayed  both  on  an  oscilloscope  or  X-Y  plotter.  A  data  flow  schematic  for  the  ground  station  is  shown  tn 
Figure  8.  Digital  tapes  for  post -flight  processing  could  be  written  simultaneously  during  data  acquisition  or  afterwards  by  playback 
of  PCM  data  from  the  analogue  tapes. 

2.4  Test  Procedure 

The  basic  procedure  adopted  was  to  excite  the  airframe  by  means  of  control  pulse  inputs  generated  by  stick  raps.  The  air¬ 
frame  response,  as  sensed  by  the  accelerometers,  was  then  monitored  and  analyzed  to  determine  frequencies  and  damping  of  the 
modes  of  interest. 

The  sequence  of  the  various  flight  tests  was  chosen  in  such  a  manner  that  the  chance  of  encountering  flutter  dunng  a  flight 
would  become  greater  as  the  test  program  proceeded.  In  this  way,  it  was  assured  that  the  test  team  would  have  acquired  a  certain 
level  of  experience  before  entering  the  more  critical  configurations  For  the  LAU-5003/A  trials,  tests  were  performed  at  an  altitude 
of  approximately  7,000  ft  above  sea  level.  To  minimize  fuel  used  in  achieving  the  required  test  conditions,  the  flight  technique 
adopted  was  to  trade-off  altitude  for  speed;  pulling  up  after  each  test  point  and  diving  to  achieve  the  next.  When  a  given  test  point  was 
stabilized,  the  pilot  was  then  given  the  instruction  to  initiate  a  lateral  stick  pulse.  After  the  stick  rap,  the  pilot  remained  hands-off  for 
a  sufficient  time  for  the  response  to  decay  before  recovering.  If  the  response  records  were  satisfactory,  the  pilot  was  cleared  to  proceed 
with  the  longitudinal  stick  pulse  followed  by  a  speed  reduction  prior  to  clearance  to  the  next  speed  increment. 

The  limitations  of  flight  duration  was  always  a  p.imary  constraint,  especially  when  obtaining  the  higher  speed  points.  The 
test  aircraft  was  at  all  times  accompanied  by  a  safety/chase  aircraft.  Since  this  aircraft  was  also  required  for  other  test  flying,  it  was  not 
equipped  with  external  tanks.  Consequently  the  higher  duration  available  when  the  test  aircraft  was  carrying  inboard  fuel  tanks  was 
not  fully  realized  and  the  flight  time  was  limited  by  chase  aircraft  fuel  limits. 

The  aircraft  was  flown  in  a  ‘dogbone’  pattern  as  shown  in  Figure  9.  With  the  out  and  back  headings  shown  in  the  figure, 
optimum  telemetry  reception  could  be  achieved  except  when  turning.  This  pattern  was  later  modified  by  extending  it  directly  over 


the  ground  station  and  out  to  a  similar  distance  to  the  southwest.  In  spite  of  a  brief  loss  of  signal  when  the  aircraft  was  over  the  ground 
station,  this  modification  resulted  in  improved  flight  efficiency. 
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2.5  Measurement  Procedure 

Representing  the  aerodynamic  forces  generated  by  the  stick  pulses  as  an  impulse  or  delta  function,  the  resulting  response  of 
the  aircraft  can  be  represented  by  a  decaying  time  function  as  shown  in  Figure  10a.  For  single  mode  response,  frequency  and  damping 
can  be  obtained  using  the  usual  logarithmic  decrement  method.  More  commonly,  however,  the  response  of  the  aircraft  is  composed  of 
several  modes  of  vibration,  and  spectral  analysis  must  be  used  instead.  Figure  10b  shows  how  damping  can  be  obtained  using  the  half- 
Dower  and  central  frequency  method. 

Based  on  analytical  predictions  and  on  previous  experience,  four  structural  modes  were  of  primary  interest  in  the  LAU-5003/A 
.flutter  trials.  They  were  the  symmetric  and  anti-symmetric  wing  bending  and  torsional  modes.  To  improve  identification  of  individual 
modes,  power  spectral  densities  were  obtained  for  four  different  linear  combinations  of  the  four  accelerometer  signals.  Each  linear 
combination  was  aimed  at  enhancing  the  power  spectrum  of  one  of  the  four  modes  of  interest,  and  the  signal  summations  used  are 
summarized  in  Figure  11.  These  summations  could  be  done  readily  on  the  HP  5451B  analyzer  since  the  machine  could  be  programmed 
to  do  arithmetic  operations  on  the  four  input  signals  before  performing  the  Fourier  transform. 

Frequency  range  up  to  10  Hz  was  desired  for  the  power  spectral  densities  since  this  range  covered  the  important  fundamental 
wing  modes  of  interest.  A  data  sampling  frequency  of  20  Hz  was  selected  giving  a  frequency  resolution  of  approximately  0.02  Hz  for  a 
block  size  of  1024.  Since  the  response  of  the  wing  usually  decayed  to  essentially  ambient  level  within  approximately  3  seconds  of  the 
stick  rap,  the  analyzer  was  programmed  to  read  data  into  a  block  size  of  128  at  50  millisecond  intervals.  Positions  beyond  80  were 
then  set  to  zero  keeping  only  the  first  4  seconds  of  data.  The  block  size  was  then  switched  to  1024.  This  resulted  in  the  data  being 
stored  in  the  First  80  positions  with  the  rest  padded  with  zeros.  In  this  way,  the  data  acquisition  time  was  kept  to  a  minimum.  The 
truncation  error  was  considered  to  be  within  acceptable  limits  based  on  tests  performed  by  NLF  (Ref.  2). 

The  power  spectral  densities  were  displayed  on  an  oscilloscope  and  on  a  X-Y  plotter.  Damping  was  determined  by  manual 
location  of  the  central  frequency  and  half  power  points.  Based  on  the  damping  thus  calculated,  a  decision  was  then  made  whether  to 
proceed  to  the  next  higher  speed  test  point.  For  qu.ck  took  analysis  and  for  monitoring  during  acceleration  to  the  next  test  point, 
the  time  signals  were  displayed  on  a  strip  chart  recorder  which  usually  could  indicate  qualitatively  the  damping  of  the  modes  under 
investigation. 

3.0  POST-FLIGHT  DATA  ANALYSIS 

Damping  values  obtained  from  power  spectral  density  plots  are  simple  and  direct  with  minimum  subjective  interpretations 
when  the  modes  are  adequately  excited  and  the  frequencies  are  sufficiently  far  apart.  However,  turbulence  and  buffet,  which  can  be 
quite  intense  at  the  high  Mach  number  tests,  degrade  the  response  signal  making  damping  measurements  difficult  and  uncertain.  Also, 
for  some  store  configurations  in  the  LAU-5003/A  trials,  insufficient  response  of  the  anti-symmetric  bending  mode  was  detected  since 
the  node  line  for  this  mode  ran  across  the  aileron  so  that  aileron  forces  did  not  produce  sufficiently  large  acceleration  in  this  mode  for 
analysis  purposes.  In  order  that  more  reliable  results  can  be  obtained,  an  interactive  computer  program  has  been  developed  which  is 
more  versatile  and  accurate  than  analogue  techniques.  This  program  has  been  used  mainly  for  post-flight  data  analysis,  but  it  can  also 
be  used  to  supplement  the  on-line  analogue  monitoring  of  subcritical  damping  in  flutter  flight  tests 

3.1  Damping  from  Impulsive  Input 

Considering  the  symmetric  and  anti-symmetric  bending  and  torsion  modes  to  be  the  important  ones  in  the  flutter  trials,  the 
acceleration  at  location  ‘A*  (Fig.  5)  for  a  given  excitation  function  can  be  written  as  follows: 

YAW  =  fAW  +  PaW  +  nA ,<‘>  U> 

where  fA(t)  is  a  deterministic  signal  determined  by  the  input  excitation  function,  pA(t)  is  a  random  signal  for  the  response  due  to 
turbulence  or  buffeting,  and  nA(t)  is  the  noise  in  the  measuring  equipment.  Equation  (1)  gives  the  general  form  of  the  response 
signal,  and  the  only  assumptions  made  for  later  analysis  purposes  are  that  fA,  pA  and  nA  are  uncorrelated,  and  pA  and  nA  are 
statistically  stationary.  Similar  expressions  can  be  written  for  the  signals  from  accelerometers  ‘B\ /C*  and  ‘D\ 

Treating  the  input  as  an  impulse  function,  and  neglecting  PA(t)  and  nA(t)  for  the  time  beir.g,  Equation  (1)  can  be 
written  as: 

YA<‘>  *  -  *SB> 

A  -otAR^  A 

+  A^fie  «n(coABt  +  *AB) 

+  A^Te""STtan(wSTt  +  *gT) 

+  A*Te  “ATtsin(uATt  +  *At) 


(2) 
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where  A  is  the  amplitude  of  the  vibratory  signal  of  one  mode,  w  it  the  circular  frequency,  t  is  the  time,  a  is  the  damping  and  ^  is  the 
phase  angle.  The  aubscript  *A’  in  Y  and  superscript  ‘A’  in  A  and  ^  denote  the  right  forward  accelerometer  ‘A’.  The  subscripts  ‘SB*, 
‘AB\  ‘ST*  and  ‘AT’  denote  the  first  symmetric  bendine,  anti-symmetric  bending,  symmetric  torsion  and  antisymmetric  torsion  modes 
respectively.  Vsrog  accelerometer  ‘A’  u  the  reference,  the  amplitudes  of  the  various  modes  at  *C'  are  the  same  as  those  at  ‘A*.  The  sym¬ 
metric  modes  are  in  phase  but  the  anti-symmetric  modes  are  180  degrees  out  of  phase.  At  location  4B\  the  amplitudes  for  the  bending 
modes  may  not  necessarily  be  the  same  as  those  at  *A\  but  the  phase  angles  are  in  phase.  For  the  torsion  modes,  not  only  are  the 
amplitudes  in  general  different  at  ‘A’  and  4B\  but  the  phase  angles  differ  by  180  degrees.  Making  use  of  these  conditions,  the  four 
accelerometer  signals  from  the  front  and  rear  of  the  tip  tanks  can  be  combined  to  give  the  following: 


(YB- YD) 


*  ♦  A® 


‘AB  AB  B 
AAT 


**n(wABt  +  *A 


Knowing  the  amplitude  ratio  A^,/A^T,  the  anti-symmetric  bending  mode  can  be  separated  as  shown  in  the  above  equation.  Similar 
expressions  can  be  obtained  for  the  anti-symmetric  torsion,  symmetric  bending  and  symmetric  torsion  modes  and  they  are  represented 
by  the  functions  F2,  F^  and  F4  respectively  in  Reference  4.  The  modal  frequency  and  damping  for  the  individual  modes  can  be 
obtained  by  forming  the  power  spectral  density  curves  for  these  F  functions  and  locating  the  central  frequency  and  half  power  points 
on  these  spectra.  As  shown  in  Figure  10b,  the  damping  is  given  as 


3.2  Damping  from  Exponentially  Decaying  Input 

For  input  excitation  to  the  aircraft  of  the  form  of  an  exponential  decaying  cosine  function,  that  is,  I  %  e~a  rcoecJ0t,  where 
I,  o'  and  cj0  are  the  applied  force,  decay  coefficient  and  circular  frequency  respectively,  the  power  spectrum  of  I  can  be  written  as 
(Ref.  5). 


1  +/J2(w-  u.)2  1  +  /32(oj  +  w0)2 


where  S  is  a  constant  and  0  »  l/a\  Referring  to  Figure  12,  let 


then  0  can  be  solved  in  terms  of  t?  and  oj0  by  the  following  equation: 


1  |  r — — - *  Vi 

fi  -  —  |(4ij-  2)  +  \J  161)2  ■  16*| + 1 


for  tj  >  0.9333.  From  Figure  13  which  shows  the  response  to  an  exponential  decaying  cosine  input,  the  damping  ratio  based  on 
Sy(wmax)  and  Sy(w^)  can  be  expressed  as 

.  [(i-n2)2 +*!!#!  -  (l-n2)2^ 

g2 - - — - - —  (8) 

n2sa  -  (fi2  + 

A  similar  expression  can  be  obtained  using  Sy(wmax)  and  Sy(a>2)»  with  all  subscripts  changed  to  ‘2*.  To  determine  wn,  the  folio wmg 
expression  is  used: 


-E+  Vfi2-  4DF 


The  terms  given  in  Equations  (8)  and  (9)  are  functions  of  t*>n,  tomax,  w2,  w0, 0,  S y(wmax),  Sy(wi)  and  Sy(w2).  The  expres¬ 
sions  arc  rather  lengthy  and  they  are  given  in  Reference  4. 
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3.3  Random  Input  and  Noiae 

In  Equation  (1),  the  effect  of  noiae  in  the  measuring  equipment  if  given  by  the  term  nA(t).  In  forming  the  F  function*  auch  as 
Equation  (3),  the  noise  can  be  included  by  simply  adding  a  term  nj(t),  n2(t)  etc.  to  the  expression*  for  F.  Since  the  response  of  the 
aircraft  has  no  correlation  with  equipment  noiae,  the  effect  of  nj(t)  etc.  on  the  response  spectra  can  be  accounted  for  by  the  addition 
by  terms  Sn^(cj),  S^w)  etc.  These  noue  spectra  are  usually  small  and  neglected  in  analyzing  the  flight  test  data. 

In  the  presence  of  turbulence  or  wing  buffeting,  the  exciting  force  has  a  random  component  as  represented  by  pA(t)  in  Equa¬ 
tion  (1).  The  mode  separation  procedure  using  expressions  like  Equation  (3)  can  still  be  used  for  the  deterministic  part  of  the  response 
signal.  Spectra  of  Fj,  F2  etc.  will  show  not  only  the  modes  as  given  by  these  equations,  but  also  other  modes  which  are  excited  by  the 
random  force.  Treating  the  power  spectrum  of  this  force  to  be  that  for  a  white  noise,  the  effect  of  a  random  component  in  the  input 
excitation  on  the  determination  of  the  damping  of  a  single  mode  of  vibration  system  can  be  accounted  for  by  taking  the  excitation 
power  spectral  density  to  be  the  sum  of  a  deterministic  and  a  random  power  spectral  density  since  it  is  assumed  i\i*i  the  two  com¬ 
ponents  are  uncorrelated.  If  Sw  is  the  power  spectral  density  of  a  white  noise  random  force,  then  Equations  (8)  and  (9)  can  still  be 
used  to  determine  damping  and  resonant  frequency  provided  that  the  ratio  SW/SQ  be  included  m  the  terms  0j,  02  inc^  D  as  given  in 
Reference  4. 


If  two  of  the  modes  are  close  together,  it  will  not  be  possible  to  separate  them  and  give  satisfactory  results  if  the  turbulence 
level  or  the  buffet  intensity  is  high.  In  Figure  14  suppose  the  random  load  that  generates  the  *0*  mode  has  a  white  noise  power  spectral 
density  Sw.  To  find  the  true  damping  and  resonant  frequency  for  the  ‘a’  mode  from  the  combined  power  spectral  density  p’jt,  the 
following  equations  can  be  used  based  on  measurements  from  the  LHS  of  the  power  spectral  density  curve: 

V* 

(10) 

and 


‘Bi  ~  V 


Bx  -  AA^ 


2A, 


G2wn 


+  G 


3*%  +  G4 


(ID 


Using  measurer.  v  from  the  RHS  of  the  combined  power  spectral  density  curve,  a  similar  expression  as  Equation  (10)  can  be  ob¬ 
tained  for  the  dampn.  \cept  that  terms  A^,  Bj  and  be  replaced  by  A2,  B2  and  C2.  Theae  terms  together  with  those  appearing 
in  Equation  (11)  are  funt;.'*'**  of  wmax,  w0,wn^,  ton^,  cu2,  SY(u>max),  SyCcO}),  Sy(w2),  Sw  and  g^.  Reference  4  gives  the 
expressions  for  these  terms. 

3.4  Interactive  Computer  Program  Desc/iption 

When  the  computer  program  »  first  called,  the  accelerometer  signals  YA,  Yg,  Yc  and  Yq  together  with  the  aileron  or 
elevator  position  are  displaced  on  the  sci^n  of  a  graphics  terminal.  On  examining  theae  signals,  the  user  chooses  the  beginning  and 
end  of  the  samples  to  be  analyzed.  The  experimental  data  is  sampled  at  a  fixed  frequency,  but  the  program  can  change  the  sampling 
frequency  to  a  lower  value  if  desired.  Also,  the  sample  length  can  be  increased  by  padding  the  data  block  with  zeros.  The  user  then  has 
the  option  of  applying  a  window  to  the  data.  If  a  rectangular  window  is  chosen,  the  range  of  data  where  the  window  applies  is  inputted 
from  the  terminal.  For  an  exponential  window,  the  data  is  being  multiplied  by  the  exponential  function  and  the  program 
requires  an  input  value  for  the  damping  7. 

The  next  step  involves  computing  the  functions  YA  ±  Yc  and  Yg  ±  YD  and  displaying  the  power  spectral  density  plots  of 

A  R 

of  these  functions.  From  these  plots,  the  movement  of  a  joystick  will  determine  the  ratios  of  AA1^  A°T  etc.  m  the  expressions 

for  the  F  functions.  If  the  value  of  a  particular  A  cannot  be  determined,  an  arbitrary  value  can  be  entered  with  the  joystick.  Also, 
values  of  A  can  be  chosen  at  ‘A’  and  ‘B’  so  as  to  prescribe  values  for  the  amplitude  ratios  to  be  used  as  first  estimates  in  evaluating 
the  F  functions.  These  will  generate  a  guide  for  the  operator  to  choose  the  next  set  of  amplitude  values  which  will  improve  on  the 
mode  separation.  The  decision  that  good  mode  separation  is  achieved  is  made  by  the  operator  from  observation  of  the  power  spectral 
density  curves.  The  user  then  chooses  the  mode  to  be  analyzed  by  specifying  the  F  function. 

From  the  selected  power  spectral  density  plot,  the  user  is  prompted  by  the  computer  to  specify  a  frequency  window  in  the 
vicinity  of  which  is  the  frequency  at  which  the  power  spectral  density  is  a  maximum.  Thu  is  achieved  by  the  movement  of  the 
joystick.  An  enlarged  display  of  the  spectrum  is  then  shown  on  the  terminal  screen.  The  operator  has  the  option  of  obtaining 
and  the  frequencies  at  the  half  power  points  from  an  automatic  computer  search.  This  procedure  is  used  if  the  spectrum  is  relatively 
noise  free.  The  second  option  is  to  fit  a  smooth  curve  and  locate  a>max  using  the  joystick.  The  computer  then  determines  the  half 
power  points  and  proceeds  to  calculate  the  damping  ratio. 

The  determination  of  damping  can  be  carried  out  by  two  methods.  The  user  can  choose  the  first  option  which  treats  the 
aerodynamic  forces  generated  by  a  stick  rap  to  be  an  impulse  function.  The  damping  can  readily  be  determined  from  the  location  of 
the  peak  of  the  power  spectral  density  curve  and  the  half  power  points.  The  second  option  treats  the  behaviour  of  ^he  aerodynamic 
forces  as  that  of  an  exponential  decaying  cosine  function.  The  user  is  then  requested  by  the  computer  to  input  the  valivt  of  q  and 
wQ,  and  then  proceeds  to  evaluate  0  from  Equation  (7).  Setting  Sj  *  2  and  using  and  Awj,  g  is  determined  from  Equation  (8). 
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Using  w 2  and  Au?2  corresponding  to  §2  "  2,  another  value  of  g  can  be  obtained.  The  damping  is  obtained  from  the  average  of  these 
two  values.  Equation  (9)  is  then  evaluated  for  u;n  which  is  u/,ed  in  Equation  (8)  to  give  an  improved  value  of  g.  An  iterative  procedure 
is  set  up  and  a  solution  is  obtained  when  g  is  within  1%  in  two  consecutive  iterations.  On  completion  of  damping  analysis  for  the 
modes  of  interest,  final  plots  for  the  F  functions  and  their  corresponding  spectra  are  displayed  at  the  option  of  the  user. 

In  the  interactive  computer  program  the  analysis  of  Section  3.3  for  random  input  and  noise  has  not  been  implemented  yet. 

A  separate  program  for  investigating  the  effect  of  turbulence  and  wing  buffeting  has  been  developed  and  it  can  be  interfaced  with  the 
interactive  program  rather  easily.  The  effect  of  measuring  equipment  noise  has  not  been  considered  since  it  is  usually  small  and  can  be 
neglected.  However,  implementing  this  into  the  computer  program  is  rather  straightforward.  Using  the  values  of  o>max,  Wj  ind  W2  for 
Si  *  §2  “  2,  Equations  (8)  and  (9)  can  be  used  to  investigate  the  effect  of  tile  presence  of  a  random  input  component  of  the  form 
given  in  Equation  (1)  on  the  damping  of  a  single  mode  vibrating  system.  For  two  closely  spaced  modes  denoted  as  the  V  and  ‘0’ 
modes,  the  procedure  is  to  use  the  interactive  program  to  give  an  estimate  for  the  resonant  frequencies  and  damping  ratios  for  these 
two  modes  by  treating  them  to  be  non-interfering.  For  known  white  noise  input,  Equation  (10)  can  be  used  to  calculate  an  average 
value  of  ga  based  on  ojj  and  cj2  for  the  V  mode.  is  then  determined  from  Equation  (11).  An  iterative  procedure  is 

required  to  solve  these  two  equations  until  ga  in  two  consecutive  iterations  is  within  the  specified  tolerance.  Using  the  values  of  gcr  and 
ton^,  the  same  equations  can  be  used  to  calculate  ‘0’  mode.  The  calculations  are  repeated  until  ga  and  g^  converge  to  within  the 
desired  tolerance,  which  is  1%  in  this  program. 

4.0  COMPUTATIONAL  METHOD  FOR  FLUTTER  CLEARANCE 

A  computer  code  has  been  developed  to  analyze  flutter  characteristics  of  aircraft  carrying  under-wing  stores.  To  calculate 
the  vibration  modes,  the  aircraft  is  separated  into  its  mam  structural  ‘branches’,  i.e.  port  wing,  starboard  wing,  front  fuselage,  etc 
The  cantilever  modes  of  these  primary  structu.al  branches  are  determined  First  to  obtain  the  so-called  branch  modes.  The  complete 
aircraft  normal  modes  are  subsequently  evaluated  in  terms  of  these  branch  mode  degrees  of  freedom  and  aircraft  rigid  body  motion. 

By  this  means  a  good  physical  description  of  the  aircraft  modes  is  obtained  with  a  minimum  number  of  degrees  of  freedom.  Aero¬ 
dynamic  generalized  forces  are  computed  using  either  aerodynamic  derivative  data  or  aerodynamic  influence  coefficients.  The  flutter 
equation  for  selected  equivalent  airspeeds  and  Mach  numbers  is  set  up  as  an  eigenvalue  problem  and  is  solved  for  modal  frequencies 
and  dampings  by  the  British  ‘p*  method. 

Figure  15  shows  a  flow  diagram  of  the  analysis  procedure.  Flutter  solutions  are  obtained  by  using  a  sequence  of  three 
computer  programs.  The  first  of  these  computes  the  branch  modes,  the  second  computes  aircraft  normal  modes  and  flutter  coeffi¬ 
cients  and  the  third  computes  flutter  solutions.  A  brief  description  of  the  individual  computer  modules  is  given  in  the  following 
sections.  More  details  on  the  computer  code  can  be  found  in  Reference  6. 

4.1  Branch  Modes  Analysis 

The  method  for  deriving  natural  inodes  and  frequencies  of  non-uniform  beams  uses  the  Holzer-Myklestad  approach.  The 
method  is  currently  restricted  to  beams  having  straight  flexural  exes  and  bending  in  only  one  plane.  Figure  16  shows  the  beam  axis 
system.  The  whole  aircraft  is  related  to  orthogonal  axes  having  their  origin  on  the  fuselage  centerline  at  the  station  vertically  above 
or  below  the  wing  elastic  axis.  Positive  directions  are  x  forward,  y  starboard  and  z  downward.  Rotations  about  these  axes  are  0, 6 
and  0.  The  elastic  properties  of  the  beam  are  concentrated  in  an  idealized  weightless  beam  lying  along  the  flexural  axis.  The  beam’s 
mass  is  represented  by  a  series  of  points,  having  both  mass  and  rotary  inertia,  whose  centres  of  gravity  do  not  necessarily  he  on  the 
flexural  axis.  These  points  are  connected  to  the  idealized  beam  at  the  series  of  stations  known  as  datum  points.  The  bending  and 
torsional  flexibilities  (1/EI  and  1/GJ)  of  the  idealized  beam  vary  linearly  with  distance  along  the  beam  between  adjacent  datum 
points. 

Pylons  perpendicular  to  the  flexural  axis  may  be  attached  at  any  datum  point.  Pylons  are  assumed  to  be  below  the  main 
beam  and  have  principal  elastic  axes  aligned  with  those  of  the  main  beam  (Fig.  17).  Pylon  stores  are  represented  by  rigid  bodies 
possessing  mass,  inertias  and  cross  inertias  about  all  axes  and  offsets  in  three  directions.  The  pylon  is  permitted  to  bend  both  in  and 
across  its  plane  and  to  twist  about  its  (vertical)  axis.  The  pylon  flexibilities  are  assumed  constant  over  the  depth  of  the  pylon. 

The  beam  and  pylon  damping  being  assumed  zero,  and  the  stiffness  linear,  this  model  can  execute  harmonic  motion  in  a 
finite  number  of  modes. 

The  idealized  model  is  considered  at  an  instant  of  maximum  deflection  in  an  oscillation  at  some  chosen  frequency  u>.  The 
displacements  and  forces  acting  at  a  section  of  the  beam  are  represented  by  the  vector  q.  Due  to  the  in-plane  rigidity  of  the  model 
the  only  significant  elements  of  the  vector  for  the  beam  are: 

q  9  Z  bending  deflection 
0  bending  slope 
6  rotation 
V  lateral  shear  force 
M  bending  moment 
T  torsion 

A  similar  q  vector  with  six  degrees  of  freedom  and  s**  forces  describes  the  pylon  store. 
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The  essence  of  the  Holzer-Myklestad  approach  is  to  relate  the  vector  qn+1  at  datum  point  n+1  to  the  vector  qn.  A  matrix 
equation  can  be  generated  which  relates  the  forces,  moments  and  displacements  (six  in  all)  at  the  beam  root  to  the  forces,  moments 
and  displacements  at  the  tip.  When  boundary  conditions  are  imposed,  a  determinant  can  be  evaluated  whose  zeros  occur  at  the  natural 
frequencies  of  the  beam.  The  method  is  programmed  to  locate  the  zeros  within  a  specified  frequency  range  and  to  evaluate  the  asso¬ 
ciated  beam  deflection  modes. 

For  the  CF-5  aircraft,  front  and  rear  fuselage  cantilever  branch  modes  have  been  computed  for  the  full  fuselage  fuel 
conf:guration.  The  empty  fuselage  configuration  is  assumed  to  have  the  same  cantilever  mode  shapes  as  the  full  fuel  case.  To  obtain 
branch  modes,  the  forward  fuselage  is  divided  into  13  stiffness  segments,  5  of  which  have  associated  mass.  The  aft  fuselage  is  repre¬ 
sented  by  6  segments,  4  with  mass.  The  fuselage  stiffness  and  mass  data  is  given  in  Reference  7.  Half-wing  cantilever  branch  modes  are 
evaluated  for  each  differ  nt  wing  store  configuration.  The  wing  is  modelled  by  19  stiffness  segments,  7  of  which  have  associated  mass. 
The  stiffness  and  mass  values  used  are  also  given  in  Reference  7. 

4.2  Calculation  of  Aircraft  Modes  and  Flutter  Coefficients 
4.2.1  Aircraft  Modes 

The  method  for  deriving  aircraft  modes  takes  as  starting  point  the  availability  of  normal  branch  modes  for  aircraft  components 
such  as  wing,  front  fuselage,  rear  fuselage.  Once  these  normal  branch  modes  are  obtained,  the  derivation  of  the  complete  aircraft 
modes  is  routine  from  the  initial  weight,  stiffness  and  geometry  data  through  to  the  required  aircraft  modes.  The  method  is  applicable 
to  all  conventional  aircraft  structures  where  such  components  as  wings,  fin,  stabilizer  and  fuselage  can  be  idealized  as  beams.  Sweep- 
back,  dihedral  and  differences  in  vertical  level  between  wings,  fuselage  and  stabilizer  can  be  taken  into  account. 

The  basic  premise  is  that  the  deformed  shape  of  the  aircraft  (in  the  modes  of  interest)  can  be  adequately  represented  by  a 
summation  over  the  original  branch  modes  and  including  certain  “ngid-body”  motions  of  the  whole  aircraft.  Limited  expenence 
indicates  that,  if  the  first  k  aircraft  modes  are  required,  it  is  necessary  to  use  a  total  of  at  least  2k  branch  modes  and  rigid  body  modes. 
To  obtain  good  accuracy  the  use  of  3k  modes  is  recommended. 

For  most  aircraft  components  the  mode  of  deformation  involves  flexible  motion  m  only  3  of  the  freedoms,  any  motion  in  the 
other  3  freedoms  being  effectively  rigid  body  motion.  For  example,  wing  bending  and  torsion  affect  only  the  heave,  roll  and  pitch 
(z,  <f>, 0)  of  any  bay  on  the  wing.  For  the  forward,  sideslip  and  yawing  motions  (x,  y,  ^ )  the  whole  wing  may  be  represented  by  a 
rigid  body.  For  the  aircraft  fuselage,  the  number  of  degrees  of  freedom  is  restricted  to  3  (or  ever  2)  by  considering  separately  the 
symmetric  and  antisymmetric  aircraft  modes. 

Having  formed  the  generalized  inertia  matrix  [a],  the  corresponding  generalized  stiffness  matrix  (e)  is  derived  by  multiplying 
the  diagonal  elements  of  [a]  by  the  corresponding  values  of  where  cj  is  the  branch  mode  frequency.  The  equation  of  motion 
(with  zero  damping  or  excitation)  can  be  written  as: 


1*1*5  +  [e]  -q  »  0  (12) 

where  q,  an  n  X  1  matrix,  represents  the  amolitude  in  each  of  the  n  modes.  For  motion  in  a  normal  mode  of  vibration  i  this  becomes 

M-q,  -  M-w,2-*,  (13) 

where  4  is  the  circular  frequency  in  the  normal  mode.  The  problem  thus  reduces  to  the  determination  of  the  eigenvalues  (latent 
roots)  ana  vectors  of  Equation  (13). 

Now,  in  general,  where  rigid  body  motions  of  the  whole  aircraft  have  been  included  in  the  n  modes  chosen,  there  will  be  a 
corresponding  number  of  zero  roots  to  the  above  equation.  The  elimination  of  these  roots  facilitates  computation  of  thoee  remaining. 

Equation  (13)  can  be  p.  rtitioned  into  the  rigid  body  and  elastic  modes  (denoted  by  o  and  e  subscripts,  respectively).  Thus 

% 

14) 

=  leoel  "  0.  Equation  (14)  can  be  sim- 

(15) 

(16) 
(17) 

Equation  (15)  gives  the  latent  roots  (1/w2)  and  vectors  (q*).  The  complete  q  matrix  can  be  obtained  from  its  partitioned  components 
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The  stiffness  associated  with  the  rigid  body  modes  is  zero,  i.e.  [eOQ]  c  [eeol 
plified  to  the  following: 

q*  -  to2  ffi  Vie«'1)q* 

where 

1*  *  qe 

•“<*  t>l  "  l«eel  -  UeJ  Kol'1  l*oel 
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4.2.2  Flutter  Coefficient*  from  Strip  Derivative* 

The  lifting  surface  under  consideration  is  divided  into  strips  o?  width  AS  a*  shown  in  Figure  18.  For  two  degrees  of  freedom 
in  heave  and  pitch,  the  lift  and  moment  for  a  strip  of  the  wing  surface  oscillating  at  a  frequency  w  can  be  written  as: 

AL  -  pV2  cAS  Jit*  +  ip«j)  •  |  +  (6a  +  ii>e„)  *  aj  <18) 

AM  -  pV2c2AS  ^(mz  +  uim*)  •  ^  +  (mQ  +  ium^)  •  a  J  (19) 

where  v  m  wc/V,  V  being  the  aircraft  velocity  •vnd  c  is  the  average  chord  of  the  strip;  fi8,  mv  ma,  C  j,  mz  and  are  the  local 

aerodynamic  strip  derivatives  and  their  values  for  the  CF-5  ‘clean  wing’  aircraft  are  given  in  Reference  7.  Using  Equations  (18)  and 
(19).  the  generalized  aerodynamic  force  can  be  obtained  in  terms  of  the  strip  derivatives  as 


where 


Qm  -  “pV2  (C  +  ii>mB)  •  qn 
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(20) 
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(22) 


Qm  is  the  generalized  aerodynamic  force  in  the  r^1  mode  due  to  a  displacement  qn  in  the  nth  mode,  and  mean  frequency 

parameter  um  *  uCm/C. 

The  generalized  matrix  equation  which  is  solved  to  obtain  flutter  roots,  frequencies  and  speeds  is 

(A  +  oA)  q  +  (B^oVg  +  D)q  +  (CVg  +  E)q  *  0  (23) 

where  o  =  p/p  ,  p0  being  the  sir  density  at  sea  level  and  VE  is  the  ratio  of  the  equivalent  airpeeed  to  a  reference  airspeed  (usually 
taken  as  1000  ft/sec).  A  is  the  structural  mass  inertia  matrix,  A  the  aerodynamic  inertia  matrix,  B  the  aerodynamic  damping  matrix 
given  by  Equation  (21),  C  the  aerodynamic  stiffness  matrix  given  by  Equation  (22),  D  the  structural  damping  matrix  and  E  the 
structural  stiffness  matrix. 

4.2.3  Flutter  Coefficient*  from  Aerodynamic  Influence  Coefficient* 

The  strip  aerodynamics  method  described  above  for  flutter  calculations  u*e*  derivative*  avaiUble  for  the  ‘clean  wing’  case. 

To  include  the  effects  of  stores,  a  more  complete  aerodynamic  model  is  required  to  give  better  result*.  The  doublet-lattice  computer 
code  H7WC  (Ref.  8)  has  been  modified  to  produce  aerodynamic  influence  coefficients  (AIC’s)  for  wing/store  configurations  (Ref.  9) 
in  calculating  flutter  coefficients. 

The  evaluation  of  the  aerodynamic  coefficient*  for  the  flutter  equation  from  aerodynamic  influence  coefficient  (A1C)  data 
can  be  determined  from  the  generalized  aerodynamic  force*.  It  i»  a**umed  that  the  AIC  matrix  i*  such  a.  to  give  the  lift  and  moment 
at  selected  points  resulting  from  unit  heave  and  pitch  accelerations  of  those  points.  For  such  an  AIC  matrix  [  F ) ,  a  typical  term  in  e 
generalized  aerodynamic  force  matrix  is 

Qi,- 2-o- tap-lF] -I*)}-*,  (24) 

-  2-a-w2  •  |g,)T  •  [F]  •  {gj}  •  qj 

where  { gj )  and  ( gj )  are  vectors  of  heave  and  pitch  displacements  in  the  ith  and  1th  modes.  [  F  ]  is  a  function  of  Mach  number  and 
reduced  frequency  K  defined  by 

K  m  (25) 

K  2V 


and  C„  is  the  AIC  reference  chord.  Using  Equations  (20)  and  (24).  the  aerodynamic  damping  and  stiffness  matrices  can  be  obtained. 
Sensitivity  of  these  matrices  to  the  K  value  for  which  the  AIC  matrix  has  been  derived  can  be  reduced  if  the  AIC  matrix  has  been 
evaluated  for  two  fairly  well  separated  values  of  K.  by  including  the  effect  of  an  aerodynamic  inertia  coefficient.  The  expression 
for  generalized  aerodynamic  force,  become* 


t 
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Q,--^-|Ctt+—  -By-—  -  A, 


If  t F] real  **  exPre8#e<I  »  *  parabolic  function  1/K 


t F I REAL  ’  2  lFd  +  (Fa1 


(fa1  •  l«jl] 


y  a  C  2 


[(Ki)T  *  tFcl  ’  Ujl] 


[  F  ]  jmag  can  “  a  line®*  function  of  1/K  as 


If1imag  '  ^  *  Ifb1 


and  By  is  then  given  by  the  expreaeion 


U  PoCBCM 


[l«i)T  ‘  tFBl  •  le,l] 


a  form  compatible  with  expressions  (26)  and  (27). 

Values  for  [Fa]  and  [Fc]  in  Equation  (27)  are  obtained  from  the  following  expressions 


. ,  '%W  • K!  ■  ifk. W  • K? 

(K|-K?) 

I'd  ■  •  Kw  ■  M 


where  [Fr  ),  [Fj^]  Me  ro®*1**  at  and  &2  respectively. 

4.3  Flutter  Solutions 

Flutter  solutions  are  obtained  from  Equation  (23).  The  basic  input  to  the  program  are  the  generalized  structural  mass  matrix, 
generalized  aerodynamic  damping  matrix,  generalized  aerodynamic  stiffness  matrix,  and  generalized  stiffness  matrix.  The  aerodynamic 
inertia  and  structural  damping  matrices  can  also  be  read  into  the  program,  but  they  are  usually  omitted  in  the  calculations.  Define 

X  “  A  +  ok 

Y  -  B  ^VE  +  D  (34) 

Z  -  CV|  +  E 


Equation  (23)  can  be  written  as 


which  is  an  eigenvalue  problem  of  the  form 


Xq  +  Yq  +  Zq  *■  0 


X  m  [-- 


i 
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where 


O 

•  1 

[T]  = 

i 

-X^Z 

-5ClYJ 

(37) 


The  eigenvalues  X  are  extracted  from  the  T  matrix  using  the  Upper  Hessenberg  method. 


5.0  SOME  RESULTS  FROM  LAU-5003/A  ROCKET  LAUNCHERS 

For  illustration  purposes,  some  results  for  the  LAU-5003/A  rocket  launchers  are  given  in  this  section.  The  configuration 
corresponds  to  that  for  Mission  1  in  the  Cold  Lake  Flutter  trials'  two  LAU-5003/A  rocket  launchers  carrying  nineteen  C14  rockets 
with  Mk  1  warheads  at  the  outboard  pylon  and  eleven  C14  rockets  with  Mk  1  warheads  at  the  inboard  pylon.  All  launchers  were 
equipped  with  nosecones. 

Figure  19  shows  the  geometry  of  the  aerodynamic  configuration  used  in  the  doublet-lattice  method.  The  tip  tanks  and 
launchers  are  represented  by  slender  bodies  with  conical  noses,  cylindrical  central  bodies  and  truncated  cones  at  the  ends.  The  wing 
root  is  considered  to  end  at  a  plane  of  symmetry  at  y/S  *  0.2.  The  steady  state  spanwise  load  and  moment  distributions  at  Mach 
number  M  *  0.8  are  shown  in  Figure  20.  The  computation  is  obtained  by  giving  the  model  a  rigid  body,  nose*up  pitch  and  setting 
the  frequency  parameter  to  zero.  All  moments  are  taken  about  the  lociJ  leading  edge,  and  in  the  tip  tank  region,  the  moment  is 
about  the  wing  leading  edge  extended  outboard.  The  unsteady  load  and  moment  distributions  for  unit  heave  and  pitch  oscillation 
about  mid -chord  of  the  wing  root  are  given  in  Figures  21  and  22  for  M  “  0.8  and  reduced  frequency  K  *  0.2. 

In  Figure  23,  the  node  lines  for  the  symmetric  and  anti-symmetnc  bending  and  torsion  modes  are  shown  for  the  full  tip 
tanks  case.  It  is  seen  that  the  node  line  for  the  anti-symmetric  bending  mode  runs  across  the  aileron  and  some  difficulties  had  been 
encountered  in  flight  tests  where  insufficient  response  of  this  mode  was  detected  when  stick  pulses  were  used  to  excite  the  aircraft. 

The  interactive  computer  program  has  been  used  primarily  to  perform  post-flight  data  analysis  for  flutter  trials.  Figure  24 
shows  a  typical  display  of  the  accelerometer  signals  at  the  four  positions  Yg,  Yq  and  Yq  indicated  in  Figure  5.  Also  shown  on 
top  of  this  figure  is  the  aileron  position.  In  this  particular  case,  the  aircraft  tip  tanks  were  full,  and  a  lateral  stick  pulse  was  initiated 
by  the  pilot  to  excite  the  anti-symmetric  modes.  The  aircraft  speed  was  437  KEAS  at  an  altitude  of  7382  ft  above  sea  level. 
Approximately  2400  points  per  channel  of  data  were  transmitted  by  telemetry  to  the  ground  station  and  recorded  on  tapes.  Since 
only  the  first  few  seconds  of  the  accelerometer  signals  after  the  stick  rap  are  of  interest,  it  is  desirable  to  suppress  the  noise  by 
applying  an  exponential  window  to  the  data.  This  is  simply  done  by  inputting  into  the  graphics  terminal  the  beginning  and  end  of  the 
range  of  data  points  for  analysis  and  the  damping  for  the  exponential  window.  Figure  25  shows  the  edited  data. 

In  Figure  26  the  function  F2  for  the  anti-sym metric  torsion  mode  and  its  power  spectral  density  plots  are  shown.  The  bottom 
curve  is  an  enlarged  plot  of  the  power  spectral  density.  The  curve  with  the  open  circles  is  a  curve  fit  of  the  original  data,  and  the  user 
of  the  computer  program  is  given  the  option  of  choosing  the  number  of  points.  The  peak  of  the  power  spectral  density  curve  is 
determined  by  the  movement  of  the  joystick  of  the  graphics  terminal  to  the  point  selected  by  the  operator.  The  horizontal  and  vertical 
lines  show  the  half  power  points  and  the  frequency  correcponding  to  the  maximum  peak  respectively.  By  moving  the  joystick 
to  the  two  half  power  points,  the  values  of  co1(  Awlt  w2,  and  Acj2  are  inputted  automatically  to  the  program  for  calculating  the 
damping. 

To  determine  the  spectra  of  the  aerodynamic  forces  generated  by  stick  raps,  a  panel  method  as  described  in  Reference  10 
may  be  used.  However,  this  method  may  not  be  suitable  in  the  flutter  flight  trials  described  here  since  the  technique  requires  the 
fluctuating  pressures  on  a  large  number  of  locations  on  the  wing  surfaces  to  be  measured.  An  approximate  estimate  can  be  obtained 
by  assuming  the  forces  to  be  proportional  to  the  control  surface  displacements  generated  by  a  stick  rap.  Also,  the  force  spectrum  is 
assumed  to  be  similar  to  the  displacement  spectrum.  Figure  27  shows  typical  time  histories  of  the  aileron  and  elevator  pulses  together 
with  their  spectra.  It  can  be  seen  that  treating  the  stick  rap  as  an  impulse  function  is  only  a  rough  approximation.  The  aileron  pulse 
damps  out  much  faster  than  the  elevator  pulse.  The  experimental  power  spectral  density  plots  are  obtained  after  passing  the  signals 
through  a  high  pass  filter.  Also  shown  in  the  figure  are  the  curves  obtained  from  Equation  (5)  using  experimentally  measured  values 
of  q  and  frequencies  <oQ  or  f0.  Figure  28  shows  that  the  frequency  fQ  for  the  aileron  and  elevator  pulses  can  be  treated  to  be  approxi¬ 
mately  constant  with  aircraft  speed  and  independent  on  whether  the  tip  tanks  are  empty  or  full.  The  values  of  i?  are  also  nearly  the 
same  'ithin  the  range  of  airspeeds  considered  except  for  the  elevator  pulse  at  airspeeds  below  470  KEAS. 

Damping  and  frequencies  for  symmetric  and  anti-symmetric  modes  are  shown  in  Figures  29  to  32  for  full  and  empty  tip 
tanks.  Theoretical  predictions  using  the  doublet-lattice  and  strip  derivatives  methods  are  compared  and  it  is  seen  that  the  bending 
frequencies  are  almost  identical  in  all  the  cases  while  the  strip  derivatives  method  gives  larger  values  for  the  torsion  frequencies.  In 
the  calculations,  no  structural  damping  has  been  included. 

The  analytical  damping  values  between  the  two  methods  are  reasonably  good  for  the  symmetric  modes  while  large  differ¬ 
ences  are  detected  for  the  anti-symmetric  modes.  Also  shown  in  the  figures  are  the  on-line  flight  test  and  post-flight  data.  In 
post-flight  analysis,  the  damping  values  are  calculated  from  Equation  (8)  using  <on  determined  from  Equation  (9).  Values  of  u>0  and  rj 
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needed  in  these  equations  are  determined  from  either  the  aileron  or  elevator  pulse.  The  frequencies  obtained  from  on-line  flight 
analysis  and  post-flight  analysis  are  very  close.  When  compared  to  theory,  the  symmetric  bending  modes  agree  very  well,  and  reason¬ 
able  agreement  is  obtained  for  the  other  cases. 

The  experimental  damping  values  from  the  two  methods  show  the  largest  discrepancies  in  the  anti-symmetric  torsion  mode. 
This  may  be  due  to  the  different  windows  applied  in  the  data  reduction  and  to  the  different  forms  of  the  excitation  power  spectral 
density  used  in  determining  the  damping  values.  Comparisons  with  calculations  show  that  analytic  computations  give  reasonable 
results  and  hence  are  useful  in  predicting  flutter  trends.  Shown  also  in  Figure  31  are  results  for  the  torsional  damping  taking  into 
consideration  interference  from  the  adjacent  bending  mode.  The  value  used  for  S^,  which  is  a  nondimensional  quantity  denoting  the 
magnitude  of  the  input  white  noise  spectra,  is  taken  to  be  0.15.  This  is  chosen  arbitrarily  in  order  to  indicate  the  effect  on  the 
damping  as  no  experiment  has  been  performed  to  determine  this  quantity.  Furthermore,  the  results  given  do  not  involve  any  iteration 
between  the  two  modes  as  described  in  Section  3.3  The  differences  in  damping  are  not  very  significant,  but  the  value  of  used  is 
considered  small  for  actual  conditions  when  wing  buffeting  occurs.  Figure  32  also  gives  a  comparison  for  damping  values  between  a 
purely  exponential  decaying  cosine  input  and  one  including  a  random  input  component  with  a  value  of  SW/SQ  equals  to  0.2.  Again 
representative  values  of  SW/SQ  in  flight  tests  are  not  available  and  the  value  of  0.2  used  is  simply  for  illustration  purposes  It  does, 
however,  indicate  that  at  high  turbulence  levels  or  buffet  intensities,  corrections  should  be  applied  m  determining  damping  values. 

6.0  CONCLUSIONS 

An  overview  of  the  Canadian  capability  in  flutter  clearance  on  the  CF-5  aircraft  carrying  underwing  stores  has  been 
described.  The  Aerospace  Engineering  Test  Establishment  at  Cold  Lake  offers  excellent  range  facilities  and  on-line  data  reduction 
for  flutter  trial  tests. 

The  use  of  stick  raps  to  excite  the  modes  of  vibrations  of  an  aircraft  in  flutter  flight  tests  is  simple  and  can  be  quite  effective. 
There  are.  however,  some  store  configurations  for  the  CF-5  aircraft  which  give  insufficient  response  of  the  anti-symmetric  bending 
mode.  This  is  due  to  the  fact  that  the  node  line  for  this  mode  runs  across  the  aileron  so  that  aileron  forces  do  not  produce 
sufficiently  large  acceleration  in  this  mode  for  analysis  purposes. 

The  computer  program  for  post-flight  flutter  analysis  developed  for  the  Canadian  Forces  operates  in  the  interactive  mode. 

It  can  also  be  used  to  supplement  on-line  analogue  monitoring  of  subcritical  damping  in  flutter  tests.  The  ability  of  the  program  to 
display  the  power  spectral  density  in  any  desired  frequency  range  makes  it  quite  useful  to  analyze  weakly  excited  modes.  Treating 
the  stick  rap  to  be  an  impulse  function  can  only  be  taken  as  a  rough  approximation.  The  damping  results  obtained  by  assuming  the 
spectra  of  the  forces  generated  by  stick  raps  to  have  the  same  form  as  that  for  an  impulse  function  may  be  quite  different  from 
those  where  the  forces  are  considered  to  be  represented  by  an  exponential  decaying  cosine  function.  The  differences  are  generally 
larger  for  the  symmetric  modes  since  the  elevator  pulses  usually  oscillate  a  few  cycles  before  decaying  to  zero  in  about  1  second. 

For  the  LAU-5003/A  trials,  the  frequencies  obtained  from  on-line  flight  analysis  and  post-flight  analysis  are  very  close.  The 
damping  values  show  the  largest  discrepancies  in  the  anti-symmetnc  torsion  mode.  This  may  be  due  to  the  different  windows  applied 
in  the  data  reduction  and  to  the  different  forms  of  the  excitation  power  spectral  density  used  in  the  calculations. 

Theoretical  predictions  using  the  doublet-lattice  and  strip  derivative  methods  give  almost  identical  bending  frequencies,  while 
the  strip  derivative  method  gives  larger  values  for  the  torsion  frequencies.  Agreement  in  damping  values  between  the  two  methods  is 
reasonably  good  for  the  symmetric  modes  while  large  differences  are  detected  for  the  anti-symmetric  modes.  Comparisons  with 
experimental  results  show  that  the  analytic  flutter  computational  method  gives  reasonable  results  and  is  a  useful  tool  in  predicting 
flutter  trends. 
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FIG.  6:  AIRBORNE  DATA  SYSTEM 


FIG.  7:  DATA  CHAIN  AT  AETE 


FIG.  8'  GROUND  TELEMETRY  AND  DATA  SYSTEM 


FIG.  9:  FLIGHT  TEST  PATTERN 
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SUMMARY 


A  correction  method  for  subsonic  potential  airforces  especially  for  the  transonic 
flow  regime  is  proposed.  The  airforces  are  corrected  by  modifying  the  theoretical  pressure 
coefficients  locally  with  measured  static  pressure  slopes..  Trends  of  transonic  airloads 
and  moments  with  reduced  frequency,  Mach-number  and  mean  static  incidence  are  given.  The 
application  of  corrected  airforces  in  flutter  calculations  is  described,  and  its  effect 
on  flutter  behaviour  is  analyzed.  It  could  be  confirmed  by  analysis  that  the  aerodynamic 
damping  of  the  most  important  low-frequency  vibration  modes  is  reduced  by  aerodynamic 
transonic  effects,  which  was  indicated  by  flight  flutter  test  results. 


1  LIST  OF  SYMBOLS 
c  ,  A  c 


C  L  =  L/q  * 
CM  =  M/q  • 

g 

k 

Mo,  '  Mloc 
NORA 

q 

s 

s 


r 

w  =  w' 

O 

<X 


+  iw" 


<*'  ♦  i  0C# 


l.e. 

t.e. 


pressure  coefficient,  pressure  jump 

sectional  steady  lift  coefficient  normal  to  elastic  axis 
lift  coefficient 

wing  bending  moment  coefficient  in  1st  wing  bending  mode 
aerodynamic  damping  coefficient 
reduced  frequency 
free-stream,  local  Mach  number 

abbreviation  for  the  organisations  NLR,  ONERA,  RAE  and  AVA 

stagnation  pressure 

surface 

semispan 

flutter  speed 

reference  speed 

elements  of  NLR  correction  matrix 

mean  incidence  of  aircraft 

amplitude  of  unsteady  angle  of  attack 

velocity  potential 

longitudinal  wash 

leading  edge 

trailing  edge 


2.  INTRODUCTION 

It  is  well  known  that  the  most  common  used  linear  *]at-plaoe  theories  for  predic¬ 
ting  unsteady  airloads  fail  in  the  transonic  flow  regime.;  The  reasons  for  this  fact  are 
phenomena  as  shocks  and  their  movement  with  periodic  deflections,  flow  separation,  visoous 
or  boundary  layer  effects,  non-linear  aspects,  appearance  of  bubbles  e^c,  At  high  angles 
of  attack  these  effects  become  more  pronounced  and  seem  to  move  beyond  the  point  to  be 
treated  any  longer  by  linear  calculation  methods. 

Therefore,  there  is  an  urgent  need  for  approximate  transonic  techniques  for  estimating 
3-d  unsteady  airforces.; 

A  list  of  trials  in  this  direction  -  far  from  being  complete  -  is  given  below.- 

E.C*;  Yates  {Ref.  1,  2)  presented  a  modif led-stnp-analysis  method  for  flutter  cal¬ 
culations.;  The  steady  state  spanwise  local  lift  distributions  and  spanwise  aero¬ 
dynamic  center  locations  ac  (normal  to  elastic  axis)  are  required  as  aerodynamic 
input  parameters  for  determining  the  oscillatory  lift-  and  moment  distributions.;  If 
measured  aerodynamic  data  are  used  instead  of  pure  theoretical  ones  it  is  assumed  that 
viscous  and  transonic  effects  are  incorporated  automatically,  at  least  in  terms  of 
nearly  steady  flow. 

The  NbR  correction  procedure  of  Bergh  and  Zwaan  (Ref.  3)  is  based  on  unsteady 
pressure  measurements  for  a  given  single  vibration  mode  and  assumes  that  the  data  infor¬ 
mation  can  be  transferred  to  all  other  modes.;  The  correction  expression  is  given  by 

ac/jC8rr  »  (ac/*  i*c;)U'0-(w'+  lw")  (*) 

The  correction  is  a  local  one.,  independent  of  vibration  modes  but  dependent  of 
reduced  frequency;  w  are  the  elements  of  tne  diagonal  and  complexfwjmatrix .  They  can 
be  evaluated  if  the  above  relation  is  solved  for  w,  where  the  subscribt  "corr"  has  to  be 
replaced  by  "exp":  w  ~  (~a test  mode*'  Both  magnitude  and  phase  of  theoretical 
loading  are  changed.;  Because  only  incompressible  results  are  reported  by  the  authors, 
the  extension  of  this  technique  to  transonic  flow  conditions  is  yet  unproved  and 
questionable . 

A  strip  theory  for  calculating  the  aero-damping  in  fundamental  bending  mode  is 
described  by  Lambourne  (Ref.  4).  The  theory  starts  from  given  steady  pressure  measure¬ 
ments  on  a  rigid  model  and  its  variation  with  incidence  up  to  and  beyond  buffet  onset. 

By  using  the  rate  of  change  of  spanwise  normal  force  distributions  with  steady 
incidence,  it  is  a  quasi-steady  strip  procedure.-  Contrary  to  Yates,  Lambourne  doesn't 
introduce  Theodorsen's  complex  circulation  function  limiting  the  application  of  his 
theory  to  small  values  of  reduced  frequency. 

A  new  approach  to  transonic  loads  was  given  by  Garner  (Ref.;  5,  6),  correcting  the 
theoretical  complex  longitudinal  wash  by  local  ratios  of  test  pressure  slopes  to  theo¬ 
retical  pressure  slopes.  The  ratios  are  assumed  to  be  independent  of  pressure  inducing 
modes: 


?-3 

Even  if  an  exact  2-dimensional  nonlinear  transonic  (small  pertubation)  method  is 
available,  one  has  to  adapt  for  3-dimensional  flows.  In  a  study  reported  by  Couston 
(Ref.'  7),  three  steady  corrections  were  used  to  account  for  viscosity,  effective  Mach- 
number  and  induced  stationary  angle  of  attack  besides  one  unsteady  correction  describing 
the  effective  motion  of  the  lifting  surface  or  the  effective  boundary  condition.  Since 
even  2-d  transonic  airloads  are  not  easy  to  be  determined,  such  an  analysis  together  with 
the  necessary  modifications  is  not  very  practicable.; 


3 .  DETERMINATION  OF  THEORETICAL  AIRLOADS 

For  predicting  unsteady  airloads,  the  doublet  lattice  method  (DLM)  according  to 
Rodden,  Giesing  and  K&lmcin  (Ref,  8)  was  chosen  because  local  pressures  are  obtained 
directly  and  not  via  pressure  series.;  In  principle,  the  same  acceleration  potential 
kernel  is  used  as  for  lifting  surface  methods  but  in  a  somewhat  different  notation, 
given  by  Landahl  (Ref.  9). 

Tne  DLM  is  a  finite  element  method  and  is  characterized  by  trapezoidal  panels  the 
vorticity  distribution  of  which  is  concentrated  in  the  quarter  chord  lines  and  whose 
downwash  control  points  are  situated  at  the  3/4  chord  lines.  The  wing  planform  of  the 
fighter  typ^  aircraft  and  its  subdivision  into  panels  is  given  in  Fig.'  1 


4 ,  CORRECTION  OF  UNSTEADY  PRESSURES  AND  GENERALIZED  AIR  FORCES 

The  wing  generalized  airforces  in  the  transonic  flow  regime  are  corrected  by  modi¬ 
fying  the  local  pressures  for  the  interesting  conditions  like  configuration  (clean  air¬ 
craft  or  wing  with  stores),  sweep  angle,  Mach-number,  steady  pre-inclination  of  aircraft, 
inducing  deflection  mode,  (reduced)  frequency.  Since  the  correction  is  punctual,  ampli¬ 
tude  and  phase  (or  real  and  imaginary  part)  of  integrated  forces  were  changed  as  well. 

Two  correction  procedures  were  investigated.; 


4  .,1  Multiplicative  Correction 


The  expression  for  the  first  version  reads  as  follows: 


(3) 


The  local  factors  (ratios)  depend  on  configuration,  on  Mach-number  and  on  steady 
state  pre-inclination,  and  are  assumed  applicable  to  all  pressure  generating  deflection 
modes.;  If  the  ratios  become  negative,  the  phases  of  the  theoretical  pressures  are  changed 
besides  the  amplitudes. 


4 . 2  Additive  Correction  (Incremental  Correction) 

The  second  applied  correction  version  can  be  written  in  the  following  complex  form: 

(* v - *  (i ( w-u, „  -  ( ‘ " '1— -  ^ ) 

The  correction  by  increments  depends  (via  the  term  c*  )  not  only  on  the  Mach-number 
and  the  steady  state  inclination  but  on  the  modes  and  on  the  reduced  frequency  too.; 


The  corrected  pressures  tend  with  k  o  and  {*'  *  <»  1  towards 

(  a*C,,/aK)exp  because  the  termsaC’  ♦  iAC£  and(‘>4/§^the0)  steady  cancel  each  other. 
Both  corrections  therefore  become  more  reliable  the  smaller  the  frequencies  are. 


5.;  evaluation  of  experimental  local  pressure  slopes 

Transonic  wind  tunnel  tests  on  a  1/14  scaled  entire  airplane  model  were  evaluated 
by  plotting  vs.  chord  upper  and  lower  side  pressures  at  given  sections  normal  to  the 
elastic  axis  and  at  chosen  inclinations  #L.  If  necessary »  curves  were  smoothed 
graphically.  Then  pressure  differences  were  produced  and  their  local  slopes  according 
to 

m  A  C  5) 

\  d *  /«*p 

For  analytical  treatment,  the  DLM  method  was  employed,  and  the  wing  planform  was 
divided  into  54  panels  (see  Fig.  1 } .  At  the  theoretical  collocation  points  the  experi¬ 
mentally  given  slopes  (Eq.  5)  were  determined  by  interpolation  and  extrapolation. 


6.  REMARKS 


Both  correction  methods  can  be  justified.  As  already  mentioned,  they  are  limited 
to  low  reduced  frequencies,  but  confirming  experimental  data  are  yet  lacking.  One 
example  for  which  the  correction  could  be  approved  experimentally  is  given  in  Fig*  2., 

The  multiplicative  and  additive  correction  was  applied  to  upper-side-pressures  of  the 
"NORA" -wing-model  (for  detailed  test  conditions,  see  figure  description  and  Ref.  11).  The 
comparison  of  the  test  data  with  theoretical  and  corrected  ones  show  that  an  improvement 
can  be  obtained  with  respect  to  the  in-phase  pressures,  though  the  peak  is  overpredicted 
in  both  cases.'  For  the  imaginary  pressure  distribution  the  incremental  coirection  is  only 
a  slight  improvement  of  the  linear  theory  while  the  multiplicative  correction  apparently 
fails.  For  this  reason  and  for  reasons  concerning  trend  studies  with  k  (not  described 
here)  ,  the  latter  approach  was  no  longer  applied.; 


7 .  AERO  -  RESULTS 


7 . 1  Steady  Flow  Pressure  Curve  Slopes  (Derivatives) 

Chordwise  experimental  and  theoretical  steady  pressure  slope  distributions  are  the 
basis  of  aforementioned  correction.,  A  typical  diagram  is  shown  in  Fig.  j.;  For  each  geo¬ 
metrical  arrangement  of  the  fighter  stores  and  pylons  and  for  each  parameter  configura¬ 
tion  (Mach-number  and  angle  of  attack) ,  a  set  of  9  similar  pictures  (according  the  9  DLM 
streamwise  strips)  has  to  be  produced,  out  of  which  local  differences  were  taken  and 
fitted  into  the  modified  DLM  programme.  Expressed  in  quasi-steady  terms  a  negative  value 
of  ( d&Cp/fix J  means  that  any  deflection  mode  will  be  excited  at  that  collocation  point 
or  nearby. 
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7 . 2  Corrected  Local  Pressure  Jumps 

st 

Pressure  plots  in  the  1  wing  bending  mode  along  selected  DLM  streamwise 
sections  are  given  in  Fig.  4 In  detail,  we  can  see  theoretical  pressures  and  correc¬ 
ted  pressures  according  to  Eq.  4. 

For  the  inner  chords,  say  3  and  4,  the  agreement  between  AC^torr  and  the  pre¬ 
dictions  is  fairly  good  while  along  the  outboard  chords  the  experimentally  corrected 
curves  show  a  similar  characteristic  but  partly  large  deviations  from  theory. 


7 . 3  Corrected  Integrated  Forces  and  Moments 

st 

Wing  lift  of  the  clean  configuration  in  1  wing  bending  mode  vs.  frequency 
parameter  k  is  shown  in  Fig.  5.' 

At  constant  reduced  frequency  k  we  can  observe  that  all  real  and  imaginary  parts 
in  the  sequence  "not  corrected,  clean  corrected  2°,  4.5°,  120,\.  This  sequence  doesn't 
change  with  varying  k,  since  the  in-phase  correction  term  is  independent  of  k  and  the 
out-of-phase  term  a  linear  function  of  k. 

The  corrected  values  therefore  show  the  same  behaviour  than  the  theoretical  ones 
when  plotted  versus  reduced  frequency. 

Diagram  6  shows  the  wing  bending  moment  associated  with  wing  bending  motion. 

The  variation  of  real  parts  and  imaginary  parts  with  reduced  frequency  is  uni¬ 
form  and  steady.  The  effect  of  the  different  mean  angles  of  attack  on  the  moment 
coefficient  is  the  same  as  observed  with  the  lifts  (Fig.  5.)., 

In  Fig.  7  the  wing  bending  moment  coefficient  in  wing  bending  mode  is  plotted  ver¬ 
sus  Mach-number.  It  decreases  with  increasing  Mach-number.  For  both  Mach-numbers  wind 
tunnel  tests  with  a  tank  on  the  inboard  wing  and  a  store  on  the  outboard  wing  were 
evaluated . 

8.  FLUTTER  ANALYSIS 


The  flutter  results  presented  in  this  paper  use  the  additive  correction.  Three- 
dimensional  corrected  unsteady  aerodynamic  forces  were  calculated  only  for  the  wing. 
Previous  investigations  showed  that  in  most  cases  the  influence  on  flutter  behaviour 
caused  by  unsteady  aerodynamic  forces  on  external  stores  and  store-wing  interference 
effects  are  very  small  (Ref.  10) *  Therefore  unsteady  aerodynamic  effects  due  to  stores 
are  not  considered  here. 

According  to  the  wind  tunnel  measurements,  calculations  were  performed  for 
Mach  0.9  and  0.95. 
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8 . 1  Representation  of  Aircraft  Structure 

A  simplified  analytical  model  of  the  advanced  fighter  type  aircraft  was  used  for 
the  investigations  presented  here.  The  45°  sweptback  wing  could  be  idealized  by  a  beam 
whereas  the  attachments  for  the  inboard  and  outboard  wing  pylons  were  represented  by 
individual  points.  It  could  be  demonstrated  by  preceding  analyses  that  the  flexibility 
of  fin,  taileron  and  fuselage  have  little  effect  on  store  flutter.  For  this  reason  these 
components  were  assumed  to  be  rigid.; 

Structural  damping  as  measured  in  ground  resonance  tests  was  not  taken  into  account 

here. 

Using  the  "Required  Damping  Method"  the  flutter  equation  was  solved  for  the 
equivalent  amount  of  structural  damping  g  necessary  for  harmonic  motion. 

9.;  FLUTTER  RESULTS 

The  following  discussion  will  concentrate  on  some  selected  results  where  the 
effects  of  corrected  aerodynamic  forces  are  considerable. 

9 . 1  Clean  aircraft 

First,  flutter  results  for  the  clean  configuration  are  shown  which  were  obtained 
by  evaluating  the  corresponding  wind-tunnel  tests  on  the  same  clean  configuration. 

In  Fig.,  8  the  first  two  mode  shapes  are  depicted  which  affect  the  flutter  behaviour * 

Fig.:  9  shows  the  considerable  influence  of  modified  airforces  on  flutter  speeds  at 
Mach  0.9  for  the  three  investigated  steady-state  incidences  *0  *  Z°f  12°.  it  is 

evident  that  the  danpings  of  the  wing  bending  mode  and  of  the  fuselage  bending  mode 
decrease,  if  correction  is  applied.  Comparing  the  two  modes  we  observe  that  at  high 
angles  of  attack  the  wing  bending  mode  shows  a  detrimental  instability  while  the 
fuselage  bending  flutter  instability  occours  at  a  higher  velocity.  The  frequency  of 
the  wing  bending  mode  is  increased  from  7.9  Hz  to  12.2  Hz. 

In  addition  to  this  symmetrical  analysis  an  antisymmetrical  analysis  for  the  same 
Mach-number ,  identical  configuration  and  the  same  mean  incidence  was  done  using  the 
equivalent  antisymmetrical  aerodynamics  and  natural  modes.  Because  the  results  were 
quite  similar  they  will  not  be  shown  here. 

9 . 2  Inboard  and  outboard  underwing  store  configuration 

The  second  example  deals  with  the  influence  of  transonic  airloads  in  flutter  analyses 
for  a  symmetrical  inboard  and  outboard  underwing  store  configuration.;  The  most  important 
resonance  mode  shapes  are  depicted  in  Fig.  10.  Fig.  11  illustrates  the  influence  of  cor¬ 
rected  aerodynamic  forces  at  Ma  0.9  and  Fig.  12  shows  the  effect  at  Ma  0.95. 

Comparing  the  V-g  plots  obtained  for  two  Mach-numbers  the  most  significant 
flutter  phenomena  will  be  discussed  now. 

In  both  graphs  the  lowest  flutterspeed  is  obtained  by  using  pure  theoretical 
aerodynamics.  There  are  drastic  damping  reductions  in  the  wing  bending  mode  at  high  angles 
of  attack.  This  effect  is  especially  pronounced  for  M  =  0,95. 
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Although  for  this  configuration  the  lowest  possible  flutter  speed  was  already  found 
employing  theoretical  aerodynamics,  the  results  may  become  important,  if  wing  bending 
reveals  as  the  flutter  critical  mode  like  it  was  in  the  clean  case.  Admittedly  this  guess 
is  not  proved  by  the  present  calculations. 


9 . 3  Outboard  store  configuration 

The  last  example  is  given  by  an  underwing  outboard  store  configuration  oscillating 
symmetrically  and  three  fundamental  resonance  mode  shapes  of  which  are  sketched  in 
Fig.  13. 


As  can  be  deduced  from  Fig.  14  and  Fig.'  15  the  correction  for  mean  angles  of  2° ,  4d 
and  6°  affects  only  the  bending  branch  to  a  great  extents  Additional  statements  about 
the  flutter  behaviour  do  not  result  from  these  calculations. 


1°.  CONCLUSION: 

It  snould  be  pointed  out  that  determination  of  pressure  gradients  with  respect  to 
mean  incidence  demand  high  accuracy  and  repeatability  of  tunnel  data.'  For  a  critical 
consideration  one  should  take  this  fact  into  account. 

In  addition  we  emphasize  again  that  the  present  correction  mechanism  is  not  yet 
proved  sufficiently  by  experiment  and  is  valid  only  for  low  reduced  frequencies. 


The  present  results  indicate  that  the  conservative  flutter  behaviour  with  pure 
theoretical  aerodynamics  changes  if  corrections  are  applied..  At  high  mean  incidences 
a  possible  instability  of  wing  bending  mode  was  found. 
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FIG.  1  WING  PLANFORM  DIVIDED  INTO  PANELS. 

BOUND  VORTEX  AT  LOCAL  1/4  CHORD  LINE,  DOWN  WASH 
COLLOCATION  POINT  AT  3/4  CHORD 


FIG.  3  SLOPE  OF  STEADY  PRESSURE  DIFFERENCE  AT  VARIOUS 
TEST  CONDITIONS. 

4S*  SWEPTBACK  WING,  SECTION  *T\ 

M  -  0.9 


FIG.  2 


UPPER  SIDE  PRESSURE  DISTRIBUTION  ON 
NORA-WING  ROTATING  ABOUT  INDICATED  AXIS. 

R  »  -P’/Ptotal,  •  *  “P,,/PT0T«. 
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FIG.  5  VARIATION  OF  LIFT  -COEFFICIENT  IN  SYMM.  BENDING  FIG.  6  VARIATION  OF  WING  BENDING  MOMENT  “COEFFICIENT 

WITH  REDUCED  FREQUENCY.  WITH  REDUCED  FREQUENCY. 

M  *  0.9,  CLEAN  CONFIGURATION  M  *  °<*’  CLEAN  CONFIGURATION,  SYMMETRICAL  CASE 


0.9  0.95  1.0 

M  — ► 


FIG.  7  SYMMETRICAL  WING  BENDING  MOMENT 
COEFFICIENT  VERSUS  MACH-NUMBER. 
WING  STORE  CONFIGURATION 
ADDITIVE  CORRECTION,  k  =  1 


WING  BOOING 


RJStlAGE  BOOING 


FIG,  $  CALCULATED  SYMMETRICAL  NORMAL  MODES  FOR 
CLEAN  AIRCRAFT 


FIG.  9 


VARIATION  OF  DAMPING  WITH  FLUTTERSPEED 
FOR  CLEAN  AIRCRAFT 
MACH  =  0.9 
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WING  BEfCING 


FIG.  10  CALCULATED  SYMMETRICAL  NORMAL  MODES  FOR 
AIRCRAFT  WITH  INBOARD  AND  OUTBOARD  STORE 
CARRIAGE 
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FIG.  1 1  VARIATION  OF  DAMPING  WITH  FLUTTERSPEED 

FOR  INBOARD  AND  OUTBOARD  STORE 
CARRIAGE 
MACH  -  0.9 
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FIG.  1 2  VARIATION  OF  DAMPING  WITH  FLUTTERSPEED 

FOR  INBOARD  AND  OUTBOARD  STORE 
CARRIAGE 
MACH  -  0.96 
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SUMARY 

This  paper  presents  the  technical  details  and  results  of  a  high-speed  wind-tunnel  test  program  of 
an  aeroelastic  cantilevered  transport  type  wing  with  two  pylon-mounted  engines.  The  tests  were  con¬ 
ducted  in  the  NASA-Langley  16-foot  Transonic  Dynamic  Tunnel  (TDT)  during  December  1981..  Flutter 
boundaries  were  determined  for  an  advanced  technology  supercritical  airfoil  and  a  conventional  airfoil 
of  identical  planforras,  mass  properties,  and  stiffness.'  The  test  parameters  included  different  values 
of  model  stiffness  and  wing  loading  at  various  angles  of  attack.  The  models  were  instrumented  at  span- 
wise  wing  stations  to  determine  bending  and  torsion  deflections  and  vertical  accelerations.'  At  two 
model  wing  stations,  pressure  transdu  ers  were  distributed  along  the  chord  to  record  static  and  un¬ 
steady  oscillatory  pressures  during  the  approach  to  and  onset  of  flutter.  This  paper  presents  the  test 
program  with  results  of  the  flutter  characteristics  and  selected  steady  and  unsteady  aerodynamic  data 
for  both  airfoils  at  different  angles  of  attack  for  various  Mach  numbers  and  dynamic  pressures.; 

INTRODUCTION 

For  more  than  20  years,  engineers  have  considered  high-speed  flutter  model  testing  to  be  necessary 
as  final  verification  before  first  flight  that  an  aircraft  is  free  from  flutter y  Why  has  the  flutter 
engineer  had  to  choose  such  an  expensive  and  difficult  method?  In  the  earlv  days  of  high-speed  air¬ 
craft  design,  the  aeroelastician  realized  that  non-compressible  aerodynamic  theory  and  flutter  model 
testing  in  a  low-speed  wind  tunnel  woulc  not  define  the  sensitivity  of  flutter  to  compressible-flow 
effects,  which  could  greatly  reduce  the  flutter  speed.*  The  complex  effect  of  transonic  aerodynamics  on 
flutter  appears  to  be  even  more  important  with  the  advent  of  the  supercritical  airfoil  designs.'  Recent 
model  tests  indicate  conflicting  results  as  to  the  severity  of  the  compressibility  effects  of  super¬ 
critical  airfoils  on  flutter  speed.  Several  of  the  more  important  factors  considered  have  been 
Reynolds  number  and  angle  of  attack,  which  can  vary  considerably  throughout  the  flight  envelope  (see 
References  1  and  2) . 

An  investigation  was  initiated  by  the  Lockheed-Georgia  Company  to  determine  a  practical  and  in* 
expensive  experimental  aeroelastic  model  program  that  would  provide  insight  into  the  effects  of 
Reynolds  number  and  angle  of  attack  on  various  airfoil  designs  regarding  flutter.;  AI30,  if  we  could 
determine  the  flutter  boundaries,  why  not  gather  additional  information  from  the  same  tests:  the  steady 
and  instead y  aerodynamic  behavior  as  flutter  was  approached  and  at  the  onset  of  flutters  Phis  in¬ 
vestigation  quickly  settled  around  the  type  of  model  design  that  would  achieve  these  objectives.  In 
the  early  day3  of  high-speed  testing,  flutter  modeling  technology  followed  similar  designs  used  for 
low-speed  testing:  a  single-spar  design  with  sectionalized  aerodynamic  sections  attached  to  it.*  This 
simple  and  economical  design  made  it  relatively  easy  to  match  the  stiffness,  mass,  and  aerodynamic  air¬ 
foils  of  a  full-scale  aircraft.  High-speed  models  were  tested  for  many  years,  using  this  sectionalized 
model  construction,  however  there  was  increasing  concern  that  the  unevenness  and  gaps  between  the 
sections  created  interference  with  the  formation  of  shock  waves  and  other  aerodynamic  characteristics. 
As  a  consequence,  smooth  skin  models  have  been  built  recently.;  However,  if  the  sectionalized  model 
design  with  its  advantages  -  ease  in  matching  mass  and  stiffness,  simpler  construction,  and  ease  of 
model  changes,  thereby  shortening  tunnel  test  time  -  could  be  proven  to  be  technically  satisfactory, 
this  design  would  still  be  desirable.  In  order  to  evaluate  the  effects  of  sectionalized  models,  an  in¬ 
vestigative  test  was  conducted  ir»  the  Lockheed-Georgia  Company  Compressible  Flow  Facility  (see 
Reference  3).  The  test  model  shown  in  Figure  1  has  airstream  chordwise  slots  which  represent  two 
typical  section  widths  of  a  sectionalized  model.  The  outboard  slot  near  the  wingtip  representeo  only 
the  inboard  side  of  a  typical  section,  since  the  aerodynamic  flow  characteristics  would  normally  be 
outward  and  beyond  the  pressure  ports.;  The  slot  depths  were  nominally  0.125  inch,  and  the  leading  and 
trailing  edges  were  cut  through  approximately  three-quarters  of  an  inch  from  their  respective  edge.* 
The  gap  widths,  approximately  0.05  inch,  were  also  scaled  proportionally  to  larger  model  gaps.  The 
aerodynamic  airfoil  was  of  an  advanced  supercritical  design.; 
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Figure  1.  Slotted  Wing  Test  Model  Used  for  Aerodynamic  Study 

Five  slot  configurations  were  tested  to  simulate  the  various  conditions  that  were  considered 
important  In  comparing  a  sectlonalized  modt'  with  a  smooth-skin  model,*  The  first  slot  configuration 
was  with  the  slots  open;  the  second  configuration  was  with  soft  foam  rubber  in  the  slots;  the  third  and 
fourth  configurations  were  with  fences  in  the  slots  which  were  1/8"  and  1/16”,  respectively,  above  the 
wing  surface  to  typify  the  steps  that  result  when  a  flexible  model  twi3ts;  the  fifth  configuration  was 
with  the  slots  filled  and  sanded  smooth  with  the  surface.  Each  of  the  five  slot  configurations  was 
tegted  at  four  Mach  numbers:'  0.67,  0.80,  0.875,  and  0.95;  and  at  five  angles  of  attack:,  -2°,  -1°,  -0°, 
+1°,  and  +2°  at  each  Mach  number. : 

The  model  was  mounted  in  the  tunnel  on  the  five-component  strain-gage  balance.  The  model  had  35 
pressure  ports  distributed  chordwise:  19  upper  surface  and  16  lower  surface  at  the  inboard  slotted 
section  centerline,  and  35  pressure  ports  -  18  upper  and  17  lower  -  at  the  outboard  section.  Data  were 
recorded  at  each  Mach  number  and  angle-of-attack  combination  for  each  model  configuration.  Direct  com¬ 
parison  plots  of  C  versus  chord  station  for  each  model  configuration  were  made  for  each  Mach  number 
and  angle  of  attach.  These  comparative  pressure  plots  were  nearly  identical  overlays  and  the  slight 
difference  is  attributable  to  run  repeatability  and/or  small  differences  in  test  Mach  number.'  The  test 
results  (see  Figures  2  and  3)  show  very  conclusively  that  sectlonalized  flutter  model  design  is  aero- 
dynamlcally  accurate  v4ien  testing  in  the  transonic  speed  range.  (See  Reference  3)< 

Once  these  tests  were  completed  and  the  results  were  satisfactory,  the  test  program  that  is 
described  in  this  paper  was  undertaken. 

TEST  PUN 


The  test  plan  was  to  conduct  an  extensive  multipurpose  aeroelastic  investigation  of  an  advanced- 
technology  supercritical  airfoil  and  a  conventional  airfoil,  both  of  a  transpor -  -type  wing  in  the  high 
subsonic  and  transonic  speed  range. 

Certain  aerodynamic  and  structural  conditions  have  been  questioned  with  regard  to  their  effect  on 
flutter.*  Probably  the  two  most  prominent  in  this  category  are  the  effects  of  Reynolds  number  and  wing 
deformation  or  angle  of  attack  on  the  aeroelastic  stability.  Now  that  the  supercritical  airfoi^  has 
made  its  debut  on  the  flutter  scene,  these  and  other  effects  that  may  increase  the  sensitivity  of  the 
supercritical  airfoil  to  flutter  over  the  more  conventional  airfoils  are  even  more  paramount.-  Another 
aerodynamic  condition,  and  one  that  should  not  be  considered  of  lesser  importance  than  the  first  two, 
is  the  unsteady  oscillatory  aerodynamic  effects  on  flutter  behavior  by  the  interaction  between  the 
steady  and  unsteady  flow  fields  through  the  periodic  motion  of  the  shock  waves. 


These  three  flutter  aerodynamic  effects  -  Reynolds  number,  deformation  and  unsteady  oscillatory 
aerodynamics  -  are  the  objectives  of  this  aeroelastic  investigation,  involving  important  comparisons 
between  the  conventional  and  supercritical  airfoil. 


Figure  2.  Slotted  Model  Prenure  Distribution  Overlap  Mach  0.95 
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Figure  3.  Slotted  Wing  Root  Loading  Aerodynamic  Coefficients 


A  typical  advanced  transport  wing  with  two  pylon-mounted  engines  was  chosen  as  the  test  model  (see 
Figure  4).'  The  model  design  and  scaling  are  shown  in  Figure  5.  The  model  was  a  cantilevered  wing 
mounted  to  the  tunnel  balance  system  with  a  half-body  root  fairing.  Two  airfoil  configurations  were 
picked:  one  design,  referred  to  in  this  paper  as  a  conventional  airfoil,  was  a  "peaky"  type  airfoil 
that  is  a  forerunner  of  today’s  advanced  supercritical  airfoil  and  has  been  in  use  for  the  past  15 
years;  and  the  other  is  an  advanced  highly  aft  loaded  airfoil  that  is  25%  thicker,  non-dimensionally 
for  the  same  cruise  speed  at  a  lift  coefficient  15*  higher  than  the  peaky  airfoil  (see  Figure  6).  Both 
airfoils  have  the  same  planform  geometry,  and  the  same  two  wing  pylon-mounted  engines.  The  wing  section 
mass  properties  were  easily  duplicated  and  checked.  The  wing  stiffnesses  were  the  same  for  both  models 
since,  with  the  use  of  the  sectionallzed  model  design,  one  spar  was  used  with  both  models.  Another 
spar  was  built  four  times  stiffer  than  the  first,  thus  giving  four  times  the  test  dynamic  pressures  at 
flutter  and  four  times  the  Reynolds  number. 


Figure  4.  Aeroelastic  Model  Wing  Planform  and 
Instrumentation  Layout 
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CONVENTIONAL  AIRFOIL 


supercritical  airfoil 


Figure  6.  Model  Airfoil  Profiles  for  Instrumented  Pressure  Sections 

The  model  angle  of  attack  was  changed  remotely  through  operation  of  the  tunnel  turntable  with 
planned  test  excursions  of  -2°,  0°,  and  +1-1 /2°A  The  wing  had  a  root  incidence  of  +4  ,  which  was  con¬ 
sidered  the  0°  flight  position. 

The  wing  wa3  sectionalized  with  12  sections  made  of  balsa  covered  with  fiberglass  and  attached  to  a 
solid  aluminum  spar  at  the  section  center  bridge  system.  (3ee  Figure  7).  The  model  spar  was  built  to 
the  jig  position  of  4°  root  to  tip.: 


fiberglass 

skin 


Figure  7.  Typical  Wing  Section  Construction 


The  wing  planform  and  instrumentation  layout  are  shown  in  Figure  4.  The  elastic  axis  inboard  of 
the  break  had  a  sweep  angle  of  29-4° t  and  outboard  of  the  break  23.75°.;  The  two  engines  were  mounted 
on  pylons  suspended  beneath  the  wing  at  38%  and  63%  of  the  wing  semispan.  The  wing  was  mass-balanced 
to  represent  approx"  lately  30%  fuel  to  obtain  the  desired  flutter  mode.'  Some  of  the  fuel  mass  was  used 
for  strengthening  the  wing  sections  to  carry  the  higher  wing  loads  during  angle  of  attack  studies. 
Figure  8  shows  the  model  mounted  in  the  tunnel  with  normal  engine  design  and  Figure  9  shows  the  model 
with  the  dummy  engines. 


Figure  8.  Aeroelostic  Model  with  Normal  Engine  Design 


Figure  9*  Aeroelostlc  Mode!  with  Dummy  Engines  Installed 


The  ground  vibration  tests  for  both  test  configurations  -  the  conventional  and  supercritical  air¬ 
foil  -  were  identical  in  modal  response  and  frequencies,  verifying  the  structural  similarity.;  Figure  10 
compares  the  frequencies  of  the  analysis  and  the  model  with  the  #4  spar  (minimun  stiffness  design). 
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SUPERCRITICAL  WING  SECTIONS 
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14.76 
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7.3 

3RD  BENDING 

16.7 

16.73 

- 

j  11  ENG 

Figure  10.  Vibration  Analysis  Frequency  Comparison 


MODEL  INSTRUMENTATION 

Since  the  test  plan  was  to  record  steady  and  unsteady  aerodynamic  pressures,  two  of  the  model 
sections  were  instrumented  with  surface  pressure  transducers  and  ports.  Wing  sections  6  and  9  were 
selected  as  the  instrumented  sections,  49.6%  and  82%  of  the  wing  span,  respectively.  The  pressure 
measuring  sensors  and  ports  were  chordwise  along  the  sectior  centerline,  and  between  the  section  center 
bridge  ribs  that  attached  the  section  to  the  spar.  Each  section  had  39  pressure  ports,  19  upper  and  19 
lower  surface,  with  one  on  the  leading  edge.  Each  section  also  had  8  pressure  transducers,  4  upper  and 
4  lower  surface,  to  be  used  as  reference  and  correlation  with  the  scani-valve  pressures. 

The  wing  spar  had  bending  and  torsion  strain  gages  located  at  the  root,  20%,  40%,  60%,  and  80% 
span;  accelerometers  were  located  at  approximately  25%  and  75%  chord  at  20%,  40%,  60%,  80%,  and  95% 
span.;  The  strain  gages  and  accelerometers  defined  the  static  and  dynamic  motion  of  the  wing.  The  wing 
was  mounted  to  the  tunnel  side  wall  through  the  balance  system  to  measure  lift,  drag,  and  other  total 
wing  behavior. 

Instrumentation  Calibration 

All  of  the  model  instrumentation  was  calibrated  prior  to  wind  tunnel  testing  and  recalibrated  and 
recorded  on  the  test  magnetic  tapes  after  installation.; 

Strain  Gages  -  The  bending  and  torsion  strain  gages  were  calibrated  by  applying  loads  and  torques 
to  each  spar  and  recording  the  resultant  outputs  in  mv/lb  in. 


3-6 


Scan 1 -Valves  Transducers  -  Eight  scani-valves  were  used  during  the  test.  The  scani-valves  were 
ganged  in  groups  of  four  -  one  set  to  measure  pressure  of  the  inboard  test  section  and  one  set  to 
measure  the  outboard  instrumented  section.  Each  scani-valve  measured  one  half  of  the  section  chord* 
All  of  the  scani-valve  transducers  were  calibrated  and  certified  in  the  NASA  calibration  laboratory  for 
several  static  pressures  with  their  outputs  recorded  in  mv/psi. 

Pressure  Ports  -  The  pressure  port  orifices  were  connected  to  the  scani-valve  transducers  using 
0.043-inch  ID  temflex  tubing.  Each  outboard  section  port  used  9- foot  tubing  and  each  inboard  section 
port  used  7- foot  tubings  Each  port  was  referenced  to  the  tunnel  static  pressure  by  a  common  source 
located  in  the  wing  root  fairing.  The  system  was  checked  for  tube  blockage  or  leaks  by  applying 
several  static  pressures  to  each  port  with  their  outputs  being  recorded. 

A  wind-off  dynamic  pressure  calibration  of  the  tubes  and  scani-valve  pressure  transducers  was  per¬ 
formed.  The  calibrator  consisted  of  a  pneumatic  cylinder  driven  by  an  electromagnetic  shaker  and  a 
tube  from  the  cylinder  to  be  fitted  over  each  pressure  port.  As  the  oscillation  frequency  of  the 
shaker  was  varied,  the  signals  from  the  scani-valve  pressure  transducer  and  the  calibrated  pressure 
transducer  in  the  cylinder  were  analyzed  using  a  frequency-response  analyzer.  The  magnitude  and  phase 
of  the  scani-valve  relative  to  the  cylinder  pressure  were  expressed  at  several  static  pressure  levels,. 
The  wind-on  dynamic  pressures  of  each  port  were  adjusted  in  phase  and  magnitude  to  agree  with  the 
pressure  transducers  along  the  chord  in  a  manner  similar  to  the  method  described  in  Reference  4. 

Pressure  Transducers  -  Eight  pressure  transducers  (Endeveco  piezoresistive,  differential  pressure 
transducers  with  thread  mounting)  were  used  for  each  instrumented  section;  four  along  the  upper  surface 
and  four  along  the  lower  surface.  Each  transducer  was  referenced  to  the  tunnel  static  pressure..  The 
pressure  transducers  were  calibrated  with  the  same  tube  calibrator  as  the  pressure  ports,'  Pressure 
response  was  recorded  at  one  static  pressure  and  at  an  oscillatory  -  -assure  at  10  Hz,  which  was  near 
the  expected  flutter  frequency. 

Accelerometers  -  All  accelerometers  were  calibrated  in  mv/g  in  the  calibration  laboratory.  Phasing 
was  checked  for  correct  installation  during  the  ground  vibration  test. 

Balance  Instrumentation  -  Sensitivity  factors  were  determined  from  laboratory  tests  and  checked 
after  model  installation  by  applying  known  loads  to  the  model. 

Model  Angle  of  Attack  -  The  model  root  angle  of  attack  was  obtained  by  visually  setting  the  model 
pitch  angle  by  marks  on  the  tunnel  wall  and  more  accurately  by  recording  the  output  from  an  angular 
accelerometer  calibrated  by  using  an  inclinometer. 

TEST  RESULTS  -  FLUTTER 

The  test  was  outlined  primarily  to  determine  the  flutter  characteristics  at  transonic  speeds 
including  affects  of  angle  of  attack  and  Reynolds  number  of  a  supercritical  airfoil  and  a  conventional 
airfoil;  the  secondary  test  objective  was  to  measure  steady  and  unsteady  aerodynamic  behavior  as 
flutter  was  approached. 

In  an  approach  to  3ave  the  expensively  instrumented  pressure  sections,  the  initial  model 

configuration  was  selected  as  the  supercritical  airfoil  with  regular  sections  replacing  the  instru¬ 

mented  sections,  and  the  less  stiff  spar  affording  the  best  opportunity  of  obtaining  flutter  at  the 
lowest  dynamic  pressure.. 

Unfortunately,  this  conservative  desire  to  approach  the  first  flutter  point  a3  safely  as  possible 
may  have  been  unwise  and  cost  some  valuable  test  time..  The  first  test  runs  were  made  at  very  low 
tunnel  dynamic  pressures  (below  40  psf)  by  increasing  q  and  Mach  number  along  a  constant  tunnel  total 

head  pressure  line.  At  these  low  q’s,  the  model  behaved  very  erratic,  and  although  appearing  to  be 

near  flutter  -  it  was  beating  in  several  modes  -  the  sweeps  were  stopped  at  0.9  Mach  number.  It  wa3 
thought  that,  by  observing  a  minimum  damping  in  the  predicted  flutter  mode  by  going  under  the  flutter 
boundary,  the  critical  Mach  number  could  be  determined  and  a  safer  approach  to  flutter  at  higher  q’s 
could  be  made..  However,  in  this  q  range  this  didn’t  seem  conclusive.  Several  questions  were  raised  at 
this  time;  one,  did  the  model  have  sufficient  aerodynamic  flow  transition  point;  and  two,  did  the  large 
engines  ( aerodynamically  blow-through)  cause  unrealistic  aerodynamic  forces  coupling  with  the  wing 
modes To  eliminate  these  questions,  a  strip  of  #46  grit,  0.10  inch  wide,  was  added  along  the  wing  at 
10%  chord  on  the  upper  and  lower  surfaces,  and  the  engines  were  replaced  with  dummy  engines  that 
represented  the  mass  properties  (except  roll  inertia)  but  with  minimun  aerodynamic  effectiveness  (see 
Figure  9).  After  the  initial  runs  it  also  seemed  apparent  that  violent  flutter  would  not  occur  and  it 
would  be  reasonably  safe  to  install  the  instrumented  sections  and  continue  the  test. 

With  these  changes,  the  tunnel  pressure  was  increased  to  represent  higher  dynamic  pressures,  and  at 
a  constant  tunnel  head  pressure,  the  q  and  Mach  number  were  increased  until  flutter  was  reached.  This 
procedure  was  repeated  for  several  higher  and  lower  tunnel  pressure  lines  until  the  flutter  boundary 
for  the  supercritical  airfoil  was  plotted  (see  Figure  11).' 

Angle  of  Attack  Tests 

Although  two  noteworthy  conditions  were  observed  during  this  portion  of  the  test  program,  r.o 
explanations  are  offered  at  tnis  time  as  to  the  probable  cause.  At  several  of  the  flutter  points  with 
the  wing  at  0  Fuselage  Reference  Line  (wing  root  incidence  +4°),  the  test  conditions  of  Mach  and  q 
were  backed  off,  the  angle  of  attack  of  the  wing  root  was  increased  to  +1-1  /2°  and  subsequently 
decreased  to  -2  For  both  angle-of-attack  settings,  the  Mach  and  q  were  increased  to  or  above  the 
previous  flutter  point,  and  in  both  incidences  the  flutter  mode  was  higher  daoped  and  flutter  was  not 
obtained  for  either  the  plus  or  the  minus  angle-of-attack  settings.:  This  same  phenomenon  occurred  when 
the  conventional  airfoil  was  tested.  It  should  be  noted  that  the  wingtip  deformed  under  load  and  did 
not  necessarily  change  the  same  as  the  root. 
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The  second  interesting  phenomenon  that  occurred  with  the  supercritical  airfoil  was  a  subflutter 
low-damped  boundary.;  It  was  observed  that,  during  each  sweep  up  a  constant  pressure  line,  a  very  low 
damped,  almost  neutrally  stable  area  was  passed  through.  The  low  damped  mode  had  the  same  frequency, 
and  by  visual  comparison  the  same  mode  shape  that  eventually  fluttered,  the  low  damped  boundary 
paralleled  the  flutter  boundary  but  at  a  considerably  lower  Mach  number.  This  low  damped  boundary  was 
reasonably  narrow  and  once  through  it  the  flutter  mode  became  heavily  damped  until  the  true  flutter 
boundary  was  approached.;  Figure  12,  an  online  response  of  amplitude  versus  frequency  plot,  shows  this 
condition  occurring  at  approximately  Mach  0.75  at  9.75  Hz.  It  is  interesting  to  note  on  Figure  12  that 
the  11  Hz  mode  does  not  indicate  the  same  subflutter,  low  damping  as  the  9.75  Hz  mode.  It  appears  that 
each  mode  may  have  its  own  flutter  characteristics  since  the  9.75  Hz  mode  is  outer  wing  bending  and  the 
11  Hz  mode  is  outer  wing  torsion. 


LEGEND 

— O-  CONVENTIONAL  airfoil 
NUMSER  A  SRAR  (|  X  El,  Gj) 

— 0-  supercritical  airfoil 

NUMKR  A  SPAR  (I  X  El,  GJ) 

— o-  supercritical  AIRFOIL 

NUMIER  3  SPAR  (4  XEI,  GJ) 

/////  AREA  OF  MINIMUM  DAMPING 
/W  SUPERCRITICAL  AIRFOIL 

- SWEEPS  OF  CONSTANT 

- PRESSURE  LINES 


FLUTTER  SO  UNO  ary 


Figure  11.  Flutter  Boundaries  for  Supercritical  Airfoil 
vs.  Conventional  Airfoil 


FREQUENCY  -  Hx 


Figure  12.  Approach  to  Flutter  Response  Amplitude  for  Two 
Modes  Showing  Subflutter  Minimum  Damping 
Prior  to  Flutter 


Unlike  the  angle-of-attack  flutter  Similarity  between  the  supercritical  and  conventional  airfoils, 
the  low  damped  boundary  observed  with  the  supercritical  airfoil  did  not  occur  with  the  conventional 
airfoil. 

The  flutter  boundary  for  the  conventional  airfoil  was  obtained  using  the  same  test  technique  used 
for  the  supercritical  airfoil. 

The  test  to  this  point  had  been  conducted  for  both  airfoils  using  the  same  spar  and  stiffnesses. 
Figure  11  shows  the  flutter  boundary  for  each  airfoil.  These  boundaries,  although  showing  different 
severity,  did  represent  the  same  flutter  mode  shape  and  frequency. 

Reynolds  Humber  Tests  -  The  following  test  objective  was  to  Investigate  the  effect  of  Reynolds 
number  on  flutter.  The  plan  called  for  testing  a  stiffer  spar  causing  the  flutter  boundary  to  occur 
at  higher  dynamic  pressures,  and  thus  higher  Reynolds  number.  The  stiffnesses  of  Spar  #3  were  in¬ 
creased  by  a  factor  of  four  uniformly  over  Spar  #4  stiffnesses  used  in  the  initial  tests.;  This  would 
give  a  Reynold 3  number  four  times  greater,  which  was  considered  the  maximun  practical  dynamic  pressure 
at  which  this  type  of  model  construction  could  safely  be  tested 

Constant  pressure  sweeps  were  made  for  increasing  pressures  in  an  attempt  to  locate  the  bottom  of 
the  flutter  boundary.  After  the  first  sweep,  it  was  observed  that  minimum  damping  occurred  at  the  same 
Mach  number  as  the  first  spar  tests.  The  sweeps  were  terminated  at  Mach  0.90  after  minimun  damping  was 
passed.  When  the  pressure  sweep  line  was  tested  v*iere  flutter  wa:  expected  but  did  not  appear,  it  was 
observed  that  the  engines  were  responding  violently  laterally  at  18  Hz  just  below  the  wing  flutter  mode 
of  20  Hz.  It  was  decided  to  make  an  even  higher  tunnel  pressure  sweep  in  an  attempt  to  make  the  wing 
flutter,  but  again  the  engine  lateral  vibrations  seemed  to  kill  the  wing  flutter  mode,  which  at  this 


3-8 


point  wa3  very  strong.'  The  dummy  engines  were  like  a  blade  and  had  considerable  aerodynamic  lateral 
lift,  so  much  in  fact  that  the  engines  were  almost  at  flutter To  eliminate  this  coupling 

the  engine  fairings  were  greatly  cutaway  minimizing  flat  acrooynamic  surface.'  Again,  a  pressure 
sweep  was  made  below  sr.ticipatea  nutter  point,  and  this  time  the  wing  flutter  mode  was  clean  and 

passed  through  a  minimun  damping  at  approximately  Mach  number  0.88  with  little  interference  from  the 
engine  mode.  The  third  pressure  line  sweep  netted  flutter  at  the  same  Mach  number  as  with  the  previous 
spar  and  at  a  q  four  times  higher,  indicating  little  or  no  Reynolds  number  effect. 

Before  the  discussion  of  flutter  is  left,  some  comment  should  be  made  about  the  flutter  charac¬ 
teristics  of  the  cantilevered  wing;’  conventional  and  supercritical  airfoils  both  exhibited  the  same 
aeroelastic  modal  response,  frequency,  and  mode  shape,  and  both  exhibited  two  flutter  modes.'  Prior  to 
ma3s-balancing  the  model,  a  preliminary  flutter  analysis  parametric  study  had  been  conducted.  A 
flutter  condition  was  desirable  that  had  an  attainable  flutter  speed  -  not  too  violent  or  conversely 
not  too  shallow  an  approach  to  flutter  and  with  little  coupling  with  otn*»r  medca,  From  this  analysis, 
a  wing  fuel  loading  of  30*  was  chosen  which  demrm2traleu  a  clean,  moderate  flutter  behavior  with  a 
flutter  frequency  of  11  Hz.  During  the  model  tests,  the  11  Hz  mode  which  was  primarily  outer  wing 
torsion  with  inner  wing  bending  became  the  predominate  mode.'  As  the  Mach  and  q  were  increased  and 
flutter  seemed  eminent  in  the  11  Hz  mode,  a  slightly  lower  frequency  mode,  9.75  Hz,  which  was  mostly  a 
third  bending-torsion  mode  (that  proved  less  desirable  for  recording  unsteady  pressure  measurements 
since  the  outer  node  line  was  near  the  outer  instrumented  section  center  chordline)  became  rapidly  less 
damped  and  in  every  case  fluttered  before  ♦‘he  higher  frequency  predicted  flutter  mode  (see  Figure  12)., 
Subsequent  to  the  test,  further  analysis  has  been  conducted  and  as  shown  in  Figure  13,y  the  two  flutter 
modes  cross  at-  nearly  the  same  wing  fuel  loading,  it  was  concluded  that  the  test  objectives  could 

be  accomplished  with  the  9.75  Hz  flutter  mode.  In  a  future  te3t,  it  will  be  interesting  to  obtain 
similar  data  with  the  11  Hz  flutter  modes  by  using  a  lower  fuel  loading.; 

TEST  RESULTS  -  STEADY  STATE  AERODYNAMICS 

Figure  14  is  a  map  of  the  test  points  that  are  presented  in  this  paper. ;  They  were  chosen  from  the 
many  steady  state  pressure  measurements  that  were  recorded  as  a  good  cross-section  of  the  test  data. 
Several  of  the  test  points  show  pressure  distributions  for  the  -2  ,  0°,  and  +1-1/2°  wing  root  angle  of 
attack  conditions.  The  flutter  boundaries  are  included  in  Figure  14  for  reference  to  location  of 
pressure  data  presented.  Each  pressure  point  on  Figure  14  gives  the  plot  number  and  root  angle  of 
attack  [examples  5  (4,12°)  and  11  (1,963°)].  The  scam-valve  steady  pressure  data  for  the  forward  50* 
chord  lower  surface  inboard  section  was  not  reduced  initially  due  to  online  computer  problems  but  will 
be  reduced  later. y 

The  steady  state  pressure  plots  are  grouped  in  Figure  15  for  all  the  conventional  airfoil  pressure 
distributions  and  in  Figure  16  for  all  the  supercritical  data.; 


Figure  13,  Flutter  Analysis  Subsequent  to  Test 


Figure  14.  Conditions  Presented  fa  Steady  State 
Aerodynamic  Pressure  Plot  Comparisons 
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Figure  16.  Supercritical  Airfoil  Steady  State  Aerodynamic  Pressure  Distribution 
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Figure  17,  Conventional  Airfoil  'Jniteady  Pressure  Distribution,  Mach.  No.  0.852, 
q  =  83.87  psf. 
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Figure  18.  Supercritical  Airfoil  Unsteody  Pressure  Distribution/  Mach.  No.  =0.80, 
q  =  73.902  psf. 
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Each  pressure  plot  gives  the  plot  number,  airfoil,  Mach  number,  dynamic  pressure,  Reynolds  number  * 
root  angle  of  attack  and  inboard  oi  outboard  test  section  angle  of  attack  with  reference  to  the  tunnel-' 

The  test  sections'  true  angles  of  attack  were  determined  from  the  bending  and  torsion  strain  gage 
data  at  each  of  five  spanwise  locations,  A  corresponding  set  of  applied  loads  was  derived  and  the  de¬ 
flections  and  twists  then  determined  from  simple  beam  theory  using  experimentally  verified  stiffness 
data . 

TEST  RESULTS  -  UNSTEADY  AERODYNAMICS 
Data  Reduction 

Four  14-channel  tape  recorders  were  used  to  store  the  pressure  and  load  information  in  real  time,; 
Typically,  each  recorder  had  one  of  the  four  sets  of  pressure  transducers,  the  accelerometers 
associated  with  them,  and  the  two  associated  scani-valves .  The  remaining  tape  channels  on  each 
recorder  were  used  for  the  other  accelerometers  and  strain  gage  information  necessary  for  the  calcula¬ 
tion  of  the  airfoil  twisty 

To  analyze  this  information,  a  HP5451C  Fourier  Analysis  System  with  a  64-channel  multiplexer  was 
used.  Twelve  channels  were  sampled  simultaneously,  with  the  system  being  triggered  on  the  first 
scani-valve  step.:  This  procedure  provided  assurance  that  the  proper  time  relation  between  sensors  was 
preserved.'  The  computer  then  contained  a  collection  of  time  histories,  each  of  which  represent  one 
second  of  information,  made  up  of  1024  points.  This  time  interval  gives  approximately  10  cycles  of 
data  at  the  flutter  frequency. 

Once  these  time  histories  were  stored  on  the  mass  storage  device,  calibration  factors  were  applied 
and  the  data  were  digitally  filtered.-  The  time  history  is  Fourier  transformed  into  the  frequency 
domain.  Then,  all  the  Fourier  components  except  in  the  frequency  band  of  interest  are  cleared.-  At 
this  point  the  operator  is  free  to  retransform  the  data  back  into  the  time  domain  or  to  leave  it  where 
it  is,  depending  on  the  type  of  presentation  to  be  used.' 

Data  presentation 

Two  different  presentations  are  used:  a  frequency  domain  and  a  time  domain.  The  frequency 
presentation  gives  phase  and  amplitude  with  respect  to  a  reference  versus  chordwise  position  (see 
Figure  17  conventional  airfoil,  and  Figure  18  supercritical  airfoil) < 

The  time  domain  presentation  i3  somewhat  more  complex,  since  the  program  uses  the  filtered  time 
history,:  At  3ome  point  in  the  time  history,  the  program  reads  the  frequency  and  amplitude,  and 

continues  to  do  so  at  predetermined  intervals  for  one  complete  cycle.  The  values  are  then  plotted  as 
vectors  for  the  appropriate  pressure  distribution.'  The  plots  are  at  15-degree  intervals  during  one 
cycle  (see  Figure  19  conventional  airfoil  and  Figure  20  supercritical  airfoil).  The  reference  chord 
(25  to  75%  C)  shown  at  each  15-degree  interval  represents  the  pitch  and  plunge  of  the  instrumented  test 
sections  at  that  time  frame,-  The  unsteady  pressure  data  were  recorded  during  dynamic  response  of  the 
wing,  excited  by  tunnel  turbulence,  as  flutter  was  approached.: 

COMMENTS 

There  are  several  results  from  this  test  where  the  authors  have  not  attempted  to  draw  any  con¬ 
clusions.-  It  is  the  objective  that  as  these  data  are  studied  and  analytical  correlations  made  that  a 
better  understanding  will  be  had  regarding  the  subflutter  low  damping  boundary  and  the  flutter  speed 
increases  observed  with  both  positive  and  negative  changes  in  angle  of  attack. 

Future  tests  investigating  various  flutter  modes,  stiffness  and  mass  distributions,  and  store 
locations  will  be  of  great  importance  in  studying  flutter  characteristics  for  the  developing  new  air¬ 
foils.  The  use  of  surface-mounted  pressure  transducers  exclusively  in  lieu  of  a  mixture  between 
transducers  and  pressure  ports  will  be  impor  .ant  in  obtaining  steady  and  unsteady  aerodynamic  pressures 
simultaneously  for  all  locations  up  to  and  including  flutter.  This  will  greatly  reduce  the  data  re¬ 
duction  errors  associated  with  tube  length,  phasing  and  attenuation  and  data  compared  at  different 
cycles. 
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Figure  19.  Conventional  Airfoil  Unsteady  Aerodynamic  Pressure  Distribution, 
Mach  No,  =0.852,  q  =  83.87  psf. 
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Figure  20 , 


Supercritical  Ai-roil  Unsteady  Aerodynamic  Pre»ur«  Distribution, 

Mach  No.  =0.80,  q- 73.902  prf. 
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SUMMARY 

Methods  of  adaptive  control  have  been  applied  to  suppress  a  potentially  violent  flut¬ 
ter  condition  of  a  half-span  model  of  a  lightweight  fighter  aircraft.  This  marked  the 
confluence  of  several  technologies  with  active  flutter  suppression,  digital  control  and 
adaptive  control  theory  the  primary  contributors.  The  control  algorithm  was  required 
to  adapt  both  to  slowly  varying  changes,  corresponding  to  changes  in  the  flight  cordition 
or  fuel  loading  and  to  rapid  changes,  corresponding  to  a  store  release  or  the  transition 
from  a  stable  to  an  unstable  flight  condition.  The  development  of  the  adaptive  control 
methods  was  followed  by  a  simulation  and  checkout  of  the  complete  system  and  a  wind  tunnel 
demonstration.  As  part  of  the  test,  a  3tore  was  released  from  the  model  wing  tip,  trans¬ 
forming  the  model  abruptly  from  a  stable  configuration  to  a  violent  flutter  condition. 

The  adaptive  algorithm  recognized  the  unstable  nature  of  the  resulting  configuration  and 
implemented  a  stabilizing  control  law  in  a  fraction  of  a  second.  The  algorithm  was 
also  shown  to  provide  system  stability  over  a  range  of  wind  tunnel  Mach  numbers  and 
dynamic  pressures.; 

INTRODUCTION 

Active  flutter  suppression  has  long  been  regarded  as  a  means  of  avoiding  flutter 
placards  of  high  performance  aircraft  when  they  are  carrying  stores.;  Research  in  this 
area  has  progressed  from  analytical  studies  to  wind  tunnel  tests  (Ref.  1-2)  and  limited 
flight  tests.  This  research  has  led  to  a  high  degree  of  confidence  in  the  methods  and  the 
hardware  used  in  this  technology.  The  work  reported  on  in  this  paper  extends  this  tech¬ 
nology  in  two  important  ways:  1)  It  uses  a  digital  control  system  rather  than  the  analog 
systems  used  in  most  previous  programs,  and  2)  it  applies  methods  of  adaptive  control  theory 
to  the  active  flutter  suppression  task. 

The  use  of  a  digital  computer  to  perform  the  control  task  is  in  line  with  the  general 
transition  of  aircraft  control  systems  from  analog  to  digital  devices.  The  flexibility 
and  complexity  available  from  the  digital  computer  make  it  an  attractive  choice  for  any 
control  application.'  Adaptive  control,  with  its  requirements  for  decision-making  and  iter¬ 
ative  parameter  estimate  calculations,  makes  digital  control  a  virtual  necessity.;  The  work 
performed  to  effect  the  digital  implementation  has  been  reported  on  previously  in  Refer¬ 
ences  3  and  4.  These  references  describe  the  computer  system  used,  the  techniques  used 
to  digitize  the  control  laws  and  presents  the  results  of  a  November  1981  wind  tunnel  entry 
that  tested  non-adaptive  digital  control  laws.  This  report  stresses  the  adaptive  aspects 
oi  the  program  with  the  digital  aspects  mentioned  only  as  required. 

The  primary  motivation  for  investigating  adaptive  control  is  that  the  multiplicity 
of  stores  that  a  modern  fighter  aircraft  carries  in  numerous  combinations  makes  it  im¬ 
possible  to  anticipate  and  correct  for  all  potential  flutter  instabilities  with  s  fixed 
control  system.  The  ultimate  requirement  of  the  adaptive  system  is  that  it  be  able 
to  accommodate  any  of  these  instabilities  with  a  minimal  foreknowledge  of  the  store  con¬ 
dition.;  A  secondary  motivation  is  that  non-adaptive  control  laws  can  suffer  from  a  lack 
of  adequate  gain  and  phase  margin  at  flight  conditions  different  from  those  at  which  th^y 
v=»re  designed,;  An  adaptive  system  should  be  able  to  enhance  these  margins  and  thereby 
expand  the  aircraft's  flight  envelope. 

Before  proceeding  further,  it  is  desirable  to  define  what  is  meant  by  adaptive  con¬ 
trol  in  the  context  of  this  paper.  One  could  think  of  there  being  a  continuum  of  levels 
of  non-adaptive  and  adaptive  control,'  The  highest,  i.e.,  the  most  demanding,  level  of 
non-adaptive  control  would  entail  gain  scheduling  wherein  the  control  law  is  designed  to 
be  a  function  of  air  data  parameters  and  the  aircraft  configuration.'  The  first  level 
of  adaptive  control  is  defined  to  be  detection  and  discrimination.  At  this  level,  the 
digital  computer  is  called  upon  to  determine,  based  on  response  information,  whether  one 
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of  several  specified  conditions  exists.;  If  a  flutter  condition  is  detected,  a  control 
law  is  engaged  to  suppress  it.: 

At  the  second  level  of  adaptive  control,  it  is  assumed  that  the  parameters  of 
the  flutter  condition  are  not  known  exactly  and  input/output  data  wire  to  be  processed 
to  provide  estimates  that  result  m  an  improved  ability  to  control  the  flutter  condition 

At  a  final  level  oi  adaptability,  the  system  would  be  aFi^d  to  not  only  identify 
the  flutter  condition,  but  also  to  construct  an  appropriate  Cv  ntrol  law  to  suppress 
it  with  no  foreknowledge  of  its  nature 

The  work  performed  m  the  program  discussed  here  limited  itself  to  the  first 
two  levels  of  adaptation.  It  is  felt  that  the  third  level  exceeas  the  capability  of 
current  off-the-shelf  computers  and  possibly  the  capability  of  adaptive  algorithms. 

The  work  described  m  this  paper  was  performed  under  contract  for  the  United 
States  Air  Force  Wright  Aeronautical  Laboratories  under  a  contract  entitled,  "Test 
Demonstration  of  Digital  Adaptive  Control  Wing/Store  Flutter."  The  work  was 
performed  by  a  Northrop/Honeywell  team  in  ,'mch  Northrop,  as  the  prime  contractor, 
contributed  its  capability  in  active  flutter  suppression  (References  1  and  2) 
while  Honeywell,  acting  as  a  subcontractor,  contributed  its  expertise  on  adaptive 
control  as  it  applies  to  flutter  instabilities  (Reference  5)  The  wind  tunnel 
test  was  conducted  dunnf  April  1982  in  the  NASA  Langley  Research  Center’s  Sixteen- 
Foot  Transonic  Dynamics  Tunnel.;  AFWAL  and  NASA  engineers  participated  actively  in 
the  test  program.; 

In  this  paper,  a  theoretical  discussion  of  the  adaptive  technique  precedes 
a  description  of  the  wind  tunnel  model  and  the  instrumentation  used  to  perform 
the  test.  Limited  results  from  simulations  that  were  carried  out  prior  to  the 
wind  tunnel  test  are  given  and  are  followed  by  a  description  of  the  wind  tunnel  test 
results 

ADAPTIVE  CONCEPTS 

Two  distinct  adaptation  algorithms  were  designed  which  correspond  to  the  first 
two  levels  of  adaptation  discussed  above.;  The  first  algorithm  is  identified  as  the 
Least  Squares  Detector  and  Discriminator  (LSDD)  and,  as  the  name  implies,  uses  a 
least  squares  fitting  procedure  to  determine  whether  the  system  response  correspond 
to  one  of  several  prespecified  flutter  modes.  If  a  correspondence  is  found,  the 
second  algorithm,  the  Maximum  Likelihood  Estimator  (MLE)  is  engaged  to  provide  a 
continuous  updating  of  the  flutter  parameters  which  are  then  used  to  provide  an  appro¬ 
priate  feedback  signal  for  flutter  suppression.; 

Least  Squares  Detector  and  Discriminator 

Figure  1  is  a  block  diagram  of  the  LSDD  concept.  The  algorithm  makes  the 
decision  as  to  whether  one  of  a  number  of  possible  flutter  modes  exists  for  the  air¬ 
craft.  Detection  is  defined  as  determining  whether  any  flutter  mode  exists  while 
discrimination  decides  which  of  the  several  modes  is  present.  The  first  step  requires 
the  construction  of  signals  that  combine  the  outputs  of  the  wing  accelerometers  in 
a  way  (see  Reference  5)  that  the  resulting  signals  approximate  the  velocity  of  a 
flutter  mode,  if  it  exists.  That  is,  the  combination  strives  to  produce  a  signal 
that  has  a  Laplace  transform  of  the  form 


1  =  ,  .  Ksoie^  mi 

u  s2  +  2  CoS  +  u>2 

where  y  is  the  sensor  output,  u  is  ’be  control  surface  input  and  w,  £ ,  K  and  $  are, 
respectively,  the  frequency,  dampi'  6 ,  gain  and  phase  of  the  flutter  mode.'  The  key 
to  the  successful  performance  of  the  adaptive  algorithm  is  in  providing  an  adequate 
estimate  of  the  four  latter  parameters  when  a  flutter  condition  exists. 

The  detector  takes  the  output  of  this  combination,  after  it  has  been  passed 
through  a  bandpass  filter,  as  well  as  information  on  the  control  surface  position 
and  fits  these  data  to  a  model  of  the  form 

yk  +  Plyk-1  +  P2yk-2  =  P3uk  +  P4uk-1  <2)' 

where  k  represents  the  discrete  time  step.  The  fit  is  performed  by  finding  the  vector 
(with  components  p^,  p^/  p^  P^)  which  minimizes 

kQ  +  N 

J=  H  (yk  *  pTzk-l)2  (3) 

k  =  k0 

where  N  is  the  length  of  the  data  segment  and 

zk-l  =  (-yk-l'  _yk-2 '  V  Uk-1)T  (4! 
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The  least  squares  estimate  of  P  that  minimizes  J  is 

p = !  X>k-i  (X>x-i  yk)  (5>' 

and  the  minimum  value  of  J  is  defined  as  J*  The  accuracy  of  the  least  squares  fit  is 
defined  to  be 

kQ  ♦  N 

J*  -  jyjD,  J0  =  2  ,(yk).2  .(6)' 

k  =  k0 

Thus,  the  output  of  the  least  squares  detector  and  discriminator  is  the  "goodness” 
of  fit,  J* ,  and  the  corresponding  parameter  vector  p.  The  value  of  J*  gives  an 
indication  of  the  dominance  of  an  oscillation  and  the  accuracy  of  the  parameter  esti¬ 
mate 


The  values  of  p  can  be  used  to  calculate  the  frequency  and  damping  of  the  flutter 
mode  .  A  flutter  mode  is  said  to  exist  when  each  of  the  following  four  criteria  is 
satisfied: 

1)  JQ  is  greater  than  a  specified  value 

2)  J*  is  less  than  a  specified  value 

3)  The  frequency  estimate  is  within  a  specified  range  and 

4)  The  damping  estimate  is  less  than  a  specified  value 

If  all  four  criteria  are  satisfied,  the  MLE  is  engaged.'  It  is  seen  that  the 
definitions  of  the  criteria  include  specifications  that  can  be  thought  of  as  para¬ 
meters  m  the  algorithm.  As  the  discussion  m  the  test  results  section  of  this 
paper  indicates,  the  specification  of  these  values  was  a  major  concern  of  the  wind 
tunnel  test.  An  additional  parameter  that  required  specif ication  was  the  value  of 
N  in  Equation  3.  The  tradeoffs  involved  are  between  speed  and  accuracy,  and  between 
irritating  false  alarms,  i.e.,  the  algorithm  identifying  a  flutter  condition  when 
there  is  none,  and  the  catastrophic  failure  of  the  algorithm  to  identify  a  flutter 
condition  when  there  is  one. 

MAXIMUM  LIKELIHOOD  ESTIMATOR 

It  is  conceptually  possible  to  use  the  p  vector  of  Equation  5  to  calculate  all 
the  parameters  necessary  to  provide  adequate  control  for  the  flutter  condition.  How¬ 
ever,  it  was  found  in  practice  that  the  quality  of  the  estimates  was  insufficient  to 
provide  reliable  control.'  This  motivated  the  search  for  an  alternative  scheme  and 
the  program  of  Reference  5  identified  and  mechanized  the  MLE  technique  for  the  flutter 
suppression  application.  Figure  2  is  a  block  diagram  of  this  concept.'  The  acceler¬ 
ometer  combination  of  the  identified  flutter  mode  is  again  used  in  combination  with 
the  control  surface  input  to  obtain  estimates  of  the  flutter  parameters.  These  par¬ 
ameters  are  then  sent  to  the  flutter  controller  which  develops  a  control  law  to  suppress 
the  flutter  condition.: 

The  analytical  basis  of  the  MLE  algorithm  starts  with  the  recursive  relation 
of  Equation  2  restated  in  first  order  form: 


n  ■ 

‘  0  1 

n  * 

p3 

Wk  *  i 

-p2  -pj 

Wk 

p4  -  Pi  p3. 

with  =  (X^)^  and  initial  conditions 


where  we  have  added  two  additional  parameters  for  the  unknown  initial  conditions,;  The 
MLE  algorithm  seeks  the  values  of  the  six  p  parameters  which  minimize  a  likelihood  func¬ 
tion  that  is  a  combination  of  the  natural  logarithm  of  the  likelihood  function  based  on 
N  samples  and  a  penalty  term  added  to  keep  the  estimates  in  proximity  of  an  a  priori 
estimate  .* 
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where  is  the  residual  (y  -  generated  by  a  Kalman  filter  of  the  form 
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and  B  is  the  convanance  matrix  of  the  residual  sequence  (assumed  to  be  constant  in 
this  application) . 

The  minimization  task  is  performed  using  the  Newton-Raphson  technique  to  find 
the  increment  in  the  p  vector  which  make  the  gradient  of  L  zero: 

V2L  4P  =  -  VL  (11) 

The  stringent  time  constraints  of  the  adaptive  flutter  suppression  task  man¬ 
dates  that  a  number  of  approximations  be  made  in  the  Equation  11  calculation.  A 
key  approximation  relates  to  the  fact  that  the  v.  *s  are  linear  functions  of  p3, 
p. ,  pc  and  p^.  If  p,  and  p2  are  thought  of  as  being  fixed,  L  is  then  quadratic 
in  the  remaining  four  parameters  and  the  solution  of  Equation  11  does  not  require 
iteration.  This  approximation  was  implemented  in  the  algorithm  by  fixing  p.  and 
P2  for  five  successive  updates  of  the  through  p,  parameters  and  then  updating 
pf  through  p.  with  and  pg  held  fixed.  The  latter  update  does  require  solving 
Equation  11  for  a  perturbation  in  the  p  values  since  the  v.  are  not  linear  values 
of  and  p2 .  K 

Control  Calculation 


Once  the  six  p  parameters  have  been  determined,  tne  first  four  are  used  to 
calculate  the  physical  parameters  w,  £,  K  and  4>  of  Equation  1.  These  parameters 
are  used  in  turn  to  calculate  the  feedback  signal  that  is  used  to  control  the 
flutter  conditon.:  The  analog  form  of  the  control  signal  is  given  by: 
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(12) 


where  is  an  algorithm  input  that  specifies  the  damping  level  prescribed  for 

the  closeu  loop  system,  dc  is  a  weighting  factor  applied  to  the  damping  estimate 
and  is  a  phase  term  used  to  correct  for  phase  lags  that  are  not  estimated  by 

the  algorithm.;  These  lags  can  result  from  the  control  surface  actuator  and  from 
delays  introduced  by  the  digital  control.; 


MODEL  AND  INSTRUMENTATION 


Figure  3  is  a  representation  of  the  wind  tunnel  model  used  for  the  demonstration. 
This  model  had  been  used  previously  in  the  programs  discussed  in  References  1  and  2 
and  these  references  describe  the  model  in  some  detail.  Therefore,  only  a  basic 
description  is  given  in  this  paper. 

The  model  is  a  30  percent  scale  half  span  model  of  the  YF-17  and  is  mounted 
on  roll  bars  which  permit  rigid  body  longitudinal  degrees  of  freedom.  The  flutter 
suppression  network  basically  entails  feeding  back  compensated  accelerometer  (sensors 
in  Figure  3)  outputs  to  leading  and/or  trailing  edge  control  surfaces.  The  adaptive 
tests  were  limited  to  using  only  the  trailing  edge  surface^  Unique  features  of 
the  system  include  a  "flutter  detector,"  an  electronic  device  that  senses  when 
large  amplitude  oscillatory  responses  occur,  and  a  flutter  stopper,  which  fires  a 
mass  inside  the  AIM-7S  (Sparrow)  missile  located  on  the  outboard  pylon.  When  the 
mass  is  deployed,  the  flutter  mechanism  of  this  configuration  is  disrupted.; 


For  the  adaptive  test,  an  additional  store,  representing  an  AIM-9S  (Side¬ 
winder)  was  designed  to  be  releasable  from  the  wing  tip  launch  rail.;  When  this 
missile  was  installed  on  the  model,  the  configuration,  which  was  designated  the 
take  off  condition,  had  a  flutter  speed  that  was  very  high.  When  the  missile  was 
ejected  by  a  remote  activating  switch,  the  resulting  configuration,  designated  the 
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downloaded  condition,  had  a  relatively  low  flutter  speedy  The  release  of  the  store  at 
a  speed  above  the  flutter  speed  of  the  downloaded  configuration  imposed  a  severe  test 
on  the  adaptive  control  algorithm's  ability  to  detect  and  control  the  resulting  flutter 
condition.. 

Figure  4  depicts  the  entire  flutter  suppression  system  as  it  was  installed  at  the 
test  site.  In  the  center  of  the  figure  is  the  control  console  which  interconnected  all 
the  components  of  the  system..  The  SEL  32/55  digital  computer  executed  the  adaptive 
control  algorithm.  This  computer  is  a  true  32  bit  machine  with  a  floating-point 
processor.  Communication  with  the  wind  tunnel  model  via  the  control  console  was  per¬ 
formed  using  12  bit  analog  to  digital  and  digital  to  analog  converters.*  Other  features 
of  the  system  include  a  hydraulic  pump  to  power  the  actuators,  a  trim  controller  which 
drives  the  horizontal  tail  and  various  monitoring  instruments,  including  a  Hewlett- 
Packard  Fourier  Analyzer. 

SIMULATION 

A  key  factor  in  the  success  of  the  program  was  the  simulations  performed  prior 
to  the  wind  tunnel  entry.  Honeywell  performed  one  simulation  that  checked  out  the 
adaptive  algorithms'  performance  on  an  aircraft  model  with  a  computer  operating  in 
the  batch  mode.*  This  assisted  in  the  programming  and  debugging  of  the  algorithm 
and  provided  initial  insight  into  the  effect  of  various  algorithm  parameters. 

A  more  complete,  real  time  simulation  was  carried  out  by  Northrop  that 
included  all  the  hardware  that  was  used  in  the  wind  tunnel  demonstration.  The  be¬ 
havior  of  the  airframe  in  the  airstream  was  accounted  for  either  by  simulating 
it  on  a  twin  processor  of  the  SEL  computer  that  contained  the  adaptive  algorithm, 
or  by  using  FM  tape  data  that  had  been  recorded  during  previous  wind  tunnel  tests. 

This  simulation  served  a  number  of  functions.  For  instance,  it  assisted 
in  debugging  the  installation  of  the  algorithm  on  the  SEL  computer.*  Perhaps  the 
most  important  use  of  the  simulation  was  to  provide  information  on  the  algorithms' 
computational  requirements.  The  wind  tunnel  test  discussed  in  Reference  3  had  shown 
that  10  milliseconds  was  a  maximum  frame  time  that  could  provide  adequate  digital 
control  of  the  6.0  Hz  flutter  mode.  The  initial  version  of  the  MLE  algorithm 
was  found  to  require  a  maximum  of  30  msecs,  of  computational  time.  A  significant 
reprogramming  effort  was  required  to  distribute  the  MLE  computations  among  a 
number  of  time  frames  and  thereby  satisfy  the  10  msec  requirement. 

The  simulation  also  proved  invaluable  in  checking  out  all  the  interconnections 
required  by  the  test  equipment  and  in  adjusting  the  scale  factors  applied  in  the 
analog  to  digital  and  digital  to  analog  conversions.; 

Figure  5  shows  strip  chart  traces  from  a  simulation  which  included  a  computer 
simulation  of  the  unstable  airframe.*  The  figure  shows  that  the  LSDD  algorithm 
quickly  identified  the  unstable  condition  and  invoked  the  correct  controller.  The 
maximum  likelihood  estimator  (MLE)  then  continuously  monitored  the  sensor  output, 
computing  new  parameter  estimates  every  five  samples.  The  latest  estimates  were 
used  to  compute  revised  gain  and  phase  parameters  for  the  controller.  In  the  fig¬ 
ure,  there  is  a  perturbation  in  the  time  histories  that  occurred  when  the  control 
loop  was  opened  intentionally <  It  is  seen  that  the  algorithm  quickly  recovered 
when  the  loop  was  closed 

TEST  RESULTS 

The  wind  tunnel  test  took  place  in  April  1982  at  the  NASA  Langley  Research 
Certer's  Transonic  Dynamics  Tunnel  and  was  comprised  of  six  days  of  wind-on  testing. 
Despite  the  limited  testing  time,  a  number  of  significant  achievements  were 
demonstrated. >  Among  these  were: 

1.  The  Least  Squares  Detection  and  Discrimination  Algorithm  was  able  to  cor¬ 
rectly  identify  an  impending  flutter  condition  in  a  timely  manner.; 

2. :  The  Maximum  Likelihood  Estimator  algorithm  was  demonstrated  to  be  capable 

of  providing  rapid  estimates  of  key  parameters  and  providing  an  effect¬ 
ive  control  of  the  flutter  condition  over  a  broad  range  of  test  condi¬ 
tions  . 

3.  The  flutter  detector  and  flutter  stopper  performed  superbly,  allowing 
for  testing  up  to  a  flutter  condition  with  minimal  risk  to  the  model. 

4. ;  In  a  culmination  of  the  test  activities,  a  Sidewinder  missile  was  re¬ 

leased  from  the  wing  tip  launcher.*  The  resulting  configuration  was 
violently  unstable.  The  adaptive  algorithm  was  able  to  identify  this 
instability  and  provide  adequate  control  in  a  fraction  of  a  second. 
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A  large  proportion  of  the  testing  time  was  devoted  to  identifying  the  com¬ 
bination  of  the  large  number  of  program  parameters  that  resulted  in  the  best  algo¬ 
rithm  performance.  Once  an  acceptable  combination  was  found,  flutter  sweeps  were 
performed  that  started  at  a  subcritical  test  condition,  then  penetrated  the  flutter 
boundary  with  the  LSDD  identifying  the  flutter  behavior  and  engaging  the  MLE .  The 
sweep  was  then  continued  to  higher  dynamic  pressure  levels  until  either  the  response 
became  too  violent  for  the  safety  of  the  model  or  the  limits  on  the  control  surface 
travel  were  reached.  Figure  6  shows  a  composite  of  the  results  of  this  testing 
It  shows  that  the  adaptive  algorithm  was  able  to  stabilize  the  flutter  condition  over 
a  wide  range  of  Mach  numbers  and  dynamic  pressures* 

Figure  7  shows  a  strip  chart  record  that  provides  some  insight  into  a  key 
aspect  of  the  algorithms,  the  generation  of  the  flutter  mode  (Figure  2)  signals. 

The  first  strip  shows  the  output  of  an  accelerometer  as  it  was  obtained  from  the 
model.  The  second  strip  is  the  same  record  after  it  has  been  passed  through  an 
analog  anti-aliasing  filter  with  a  50  Hz  cutoff  frequency.  The  third  strip  is 
this  signal  after  it  has  been  digitized  while  the  fourth  trace  is  the  flutter  mode, 
obtained  by  blending  the  outputs  of  the  four  accelerometers  and  passing  the  re¬ 
sultant  signal  through  a  band  pass  filter*  It  is  seen  that  whereas  the  six  Hertz 
flutter  mode  is  masked  by  other  responses  in  the  original  accelerometer  output,  the 
final  flutter  mode  signal  is  very  clear 

The  fifth  strip  on  Figure  7  shows  the  trailing  edge  position  record  as  it 
was  measured  from  a  control  surface  potentiometer  and  the  sixth  is  the  same  signal 
after  it  has  been  digitized  and  filtered.  This  last  trace  is  the  one  used  in  combina¬ 
tion  with  the  fourth  trace  for  parameter  estimation.. 

Smoothing 

By  its  nature,  adaptive  control  makes  quantization  of  the  algroithm's  per 
formance  difficult.  Measurement  techniques  used  in  previous  studies  that  provided 
transfer  function,  power  spectral  density  and  damping  trend  data  are  not  applicable 
to  adaptive  control  testing  since  the  feedback  parameters  are  changing  constantly. 

What  follows  then  is  a  more  qualitative  assessment  of  the  test  results  than  is  nor¬ 
mally  given  in  a  paper  of  this  type. 

One  of  the  lessons  learned  during  the  tunnel  test  was  that  there  were  sig¬ 
nificant  variations  in  the  parameter  estimates  obtained  every  50  msecs  by  the 
MLE  algorithm.  These  variations  were  attributed  to  the  fact  that  the  turbulent 
environment  of  the  tunnel  gave  a  low  signal  to  noise  ratio.;  The  average  values  of 
the  estimates  appeared  acceptable  and  these  average  values  were  obtained  by  apply¬ 
ing  a  smoothing  filter  to  the  estimates.  The  smoothing  filter  used  was  a  standard 
first  order  low  pass  filter  developed  for  digital  computers  and  has  the  equation: 

Pk+i  ■  8  h  +  (1‘B)  pk+i 

where  the  k  subscript  designates  the  time  point,  the  bar  superscript  refers  to  the 
smoothed  estimate  and  6  is  the  smoothing  parameter,  which  can  vary  from  0  {no 
smoothing)  to  1  (no  updating),* 

The  selection  of  the  appropriate  B  values  became  a  major  tradeoff  activity  of 
the  wind  tunnel  test.  No  smoothing  filters  were  applied  in  the  Least  Squares 
portion  of  the  algorithm  while  five  filters  were  ultimately  applied  in  the  Maximum 
Likelihood  portion.  These  filters  can  be  classified  into  two  sets*  The  first  set 
was  incorporated  into  the  parameter  estimation  process  (Equations  7  through  11) 
and  the  second  set  was  incorporated  into  the  updating  of  the  control  law  (Equa¬ 
tion  12) 

Nominal  values  of  the  five  B 1 s  were 


Parameter (s) 

S  Value 

pr,  P2 

0.6 

p3'  ?4 

0.4 

K 

0.94 

5 

0.98 

$ 

0.92 

As  the  list  implies,  it  was  necessary  to  smooth  the  damping  estimate  by  the  largest 
amount*  The  variability  of  the  damping  estimate  was  one  of  the  more  unpleasant 
surprises  of  the  test  and  prompted  the  majority  of  the  improvisations  that  took  place 
during  the  test. 


Parameter  Variations 
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The  smoothing  filters  required  a  change  in  the  MLE  algorithm  during  the  test* 
Other  changes  were  made  during  the  test  in  the  parameters  of  both  adaptive  algo¬ 
rithms.'  A  number  of  these  parameters  will  be  discussed  briefly  here  in  order  to  give 
more  insight  into  the  algorithms. 

The  value  of  the  N  parameter  of  Equation  3  and  Equation  9  represented  a 
tradeoff  between  the  quality  of  the  estimates  and  their  timeliness.  Values  of  60 
to  100  were  used  which, with  a  sample  time  of  .01  secs,  means  that  the  data  window 
varied  from  0.6  to  1.0  sees.; 

A  parameter  denoted  SIGUT  controlled  the  amplitude  of  the  test  signal  that  was 
applied  to  the  control  surface.  In  this  case,  the  tradeoff  was  between  obtaining  a 
good  response  signal  with  a  large  value  of  SIGUT  and  keeping  SIGUT  small  enough  so 
that  the  test  signal  did  not  drive  the  model  too  hard  or  exceed  the  control  surface's 
capability.  The  parameter  was  set  so  that  the  peak  magnitude  of  the  excitation  was 
approximately  1.5°  for  the  LSDD  algorithm  and  3.0°  for  the  MLE  algorithm. 

Parameters  ZCMD  and  DC  correspond  to  and  dc  in  Equation  12.'  The  nominal 

value  of  ZCMD  was  set  at  0,02,  implying  a  request  for  2%  damping  of  the  cntical  mode 
for  the  closed  loop  system.;  The  DC  value  was  decreased  from  1.5  to  1.0  during 
the  course  of  the  test.; 

Other  parameters  controlled  the  detection  parameters  discussed  in  the  des¬ 
cription  of  the  LSDD  algorithm,  limited  the  amount  of  parameter  estimates  could 
vary,  emphasized  initial  values  of  the  parameters  to  be  estimated  and  provided 
initial  gain  and  phase  estimates  to  the  MLE  algorithm.; 

This  list  of  parameters,  which  is  still  only  a  fraction  of  those  required  by 
the  algorithm,  gives  a  sense  of  the  complexity  of  the  adaptive  algorithm.  Since  the 
wind  tunnel  test  was  the  first  opportunity  to  test  the  algorithm  in  "real  world" 
conditions,  it  is  not  surprising  that  a  large  amount  of  trial  and  error  was  required 
to  obtain  satisfactory  performance. 

Store  Release 

The  most  demanding  goal  of  the  test  was  to  demonstrate  the  ability  of  the 
adaptive  algorithm  to  detect  and  control  a  flutter  condition  that  was  entered  into 
abruptly  due  to  a  store  release.  This  goal  was  achieved  during  the  test  entry  when 
a  Sidewinder  missile  was  released  from  the  wing  tip  of  the  model  with  the  resulting 
configuration  violently  unstable.'  The  test  point  corresponding  to  this  event  is 
marked  by  a  square  in  Figure  6.'  The  LSDD  algorithm  recognized  this  instability 
very  quickly  and  engaged  the  MLE  algorithm,  which  stabilized  the  condition.'  Figures 
8  through  10  are  strip  chart  records  of  this  event  and  provide  insight  into  the  per¬ 
formance  of  the  algorithm.  Figures  8  through  10  were  reconstructed  from  the  data  that 
were  tape  recorded  during  the  test  since  the  quality  of  the  oscillograph  recordings 
obtained  during  the  test  was  poor.  The  nature  of  the  responses  in  the  reconstruc¬ 
tion  is  very  similar  to  those  obtained  during  the  test..  Figure  8  shows  variables 
from  the  digital  algorithm.;  Before  the  store  drop,  frequency  and  damping  estimates 
were  made  every  0.3  seconds,'  It  is  impossible  to  pinpoint  the  exact  moment  of  the 
store  drop  on  this  figure,  but  from  studying  all  the  data  one  can  deduce  that  there 
was  an  approximately  0.2  second  delay  between  the  store's  release  and  the  engaging 
of  the  trailing  edge  surface  control.  This  implies  that  the  LSD  algorithm  detected 
the  flutter  condition  at  its  first  opportunity  after  the  store  drop.'  The  damping 
estimate  is  seen  in  Figure  8  to  go  negative  after  which  the  MLE  algorithm  nominally 
estimates  the  open  loop  damping  to  be  a  negative  one-half  percent  of  the  critical 
damping  level.  The  large  discrepancy  between  the  magnitude  of  the  damping  estimates 
given  by  the  MLE  and  LSD  algorithms  remains  a  puzzling  factor.  Other  strips  of 
Figure  8  show  that  the  unfiltered  estimates  of  gain  and  phase  fluctuate  greatly, 
motivating  the  need  for  the  smoothing  filters  discussed  earlier.'  The  filtered 
estimates  are  seen  in  Figure  9  to  be  well-behaved.  It  should  be  pointed  out  that 
Figures  8  and  9  were  obtained  on  two  different  passes  of  the  tape  data  through  the 
adaptive  algorithm.  It  is  seen  that  the  common  strips  on  the  two  figures  do  not  cor¬ 
respond  exactly..  It  is,  however,  reassuring  that  there  is  strong  resemblance  m  the 
records,  indicating  that  the  algorithm  outputs  are  repeatable. 
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concluding  remarks 

The  test  demonstration  of  adaptive  control  of  flutter  represented  an  ambitious 
advance  in  active  control  technology.  To  the  best  of  our  knowledge ,  this  is  the 
first  time  such  a  demonstration  has  been  made  on  an  elastically  unstable  system. 

It  is  very  gratifying,  therefore,  that  the  major  goals  of  the  project,  i.e.,  the 
demonstration  of  the  LSDD  and  MLE  algorithms  and  of  the  adaptation  to  the  store 
release,  were  met  during  the  test.  This  sense  of  accomplishment  is  reinforced 
when  one  considers  that  the  test  was  carried  out  under  two  severe  time  constraints: 

1)  the  contraint  that  the  adaptive  calculations  to  be  made  within  the  10  millisecond 
frame  time  with  a  digital  computer  of  moderate  speed  performance,  and  2)  the 
constraint  that  the  wind  tunnel  test  was  limited  to  six  days. 

On  the  other  hand,  it  must  be  stressed  that  this  first  demonstration  should 
not  be  thought  of  as  the  final  word.  The  performance  rating  of  the  adaptive  te3t 
was  considered  moderate  in  both  the  quality  of  the  parameter  estimates  obtained  by 
the  algorithms  and  in  the  range  of  test  conditions  over  which  the  algorithms  stabil¬ 
ized  the  models  Improvement  in  each  of  these  areas  is  required  before  flight  tests 
using  the  adaptive  control  algorithms  can  be  safely  attempted.  The  data  recorded 
on  FM  tape  during  the  test  is  a  valuable  resource  that  can  be  utilized  in  the  search 
for  this  improvement. 

An  important  factor  in  obtaining  these  improvements  is  the  continued  rapid  advance 
in  the  capabilities  of  digital  computers.  The  most  recent  laboratory  minicomputers 
have  computational  speed  capabilities  an  order  of  magnitude  faster  than  the  computer 
used  in  this  demonstration.  Even  further  efficiencies  could  be  obtained  by  using  an 
array  processor  to  speed  the  vector  operations.  The  use  of  these  state  of  the  art 
devices  would  allow  more  detailed  and  exact  calculations  in  the  MLE  algorithm, 
with  a  resulting  improvement  in  the  parameter  estimates  and  hence  the  controller's 
capability. 

» 

While  the  computer  gains  will  provide  a  dramatic  advance  in  the  technology,  its 
viability  can  come  only  after  a  number  of  incremental  improvements  in  the  algorithm 
and  greater  understanding  of  the  basic  phenomena  with  a  flight  test  providing  the 
ultimate  demonstration. 
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Figure  l.:  Least  Squares  Detector  and  Discriminator  Adaptive 
Control  Concept. 
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Figure  4-  Flutter  Suppression  System  as  Installed  at  the  Test  Site. 
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Figure  5.  lime  History  showing  the  Adaptive 

Algorithms  Performance  When  Applied 
4-n  a  s-imuiation  of  the  Unstable 
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Figure  6.  Performance  of  the  Adaptive 
Algorithms  in  Expanding  the 
Flight  Envelope. 
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Figure  7.  Generation  of  Signals  Used  in  the  Estimation 
Process 
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Figure  8.:  Time  Histories  of  Selected  Digital  Outputs 
During  the  Score  Release. 
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This  publication  includesTeur  papers  given  to  the  Aeroelasticity  Sub-Committee  of  the 
Structures  and  Materials  Panel;  they  show  a  cross  section  of  recent  activities  by  NATO 
countries  within  the  aeroelastic  field. 

Topics  covered  are  (J)  Subsonic  flutter  clearance  procedures  in  Canada 
®)  Transonic  flutter  research  in  Germany  , 

(GJ)  Transonic  unsteady  aerodynamic  measurements  in  the 
United  States  y  >,  A 

The  use  of  adaptive  control  methods  to  actively  suppress 
flutter  of  an  aeroelastic  model  wing  with  a  store.  , 


The  papers  reflect  the  strides  being  made  in  aeroelasticity  by  the  NATO  countries. 
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